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Abstract—The mechanisms of condensation from a mixture of a vapour and a non-condensable gas 
are reviewed. A theory of condensation by mixing with a cool gas is given when the Lewis number 
for the condensing species is unity. The method can be used to predict whether fog formation is 
possible in a mixing system. An experimental system designed to test some of the theoretical 
conclusions is described, The apparatus permits the production of air-vapour mixtures of known 
and controllable temperature and concentration. Experimental data are reported for a turbulent 
jet containing glycerine vapour. The results indicate that vapour concentrations considerably 
higher than the theoretical minimum are necessary for detectable condensation in jet mixing. 


Résumé— Revue du mécanisme de la condensation d'un mélange de vapeurs et de gas inertes. 
Théorie de la condensation par mélange dans le cas de la valeur unité pour le nombre de Lewis 
de la vapeur condensable ; cette méthode permet de prédire la formation éventuelle d'un brouillard 
dans la zone de mélange. Description dun dispositif expérimental pour vérifier certaines des 
conclusions théoriques. Cet appareil permet la production et le controle de mélanges air-vapeur 
de concentrations et de températures connues., Données expérimentales pour un jet turbulent 
contenant de la vapeur de glycérine. Résultats montrant qu'il est nécessaire de dépasser con- 
sidérablement la concentration minimum théorique pour pouvoir obtenir une condensation 
décelable dans le mélange par éjecteur. 


Zusammenfassung Die Mechanismen der Kondensation eines Gemisches aus einem Dampf 
und einem nichtkondensierbaren Gas wird betrachtet. Kine Theorie der Kondensation durch 
Zumischen cines kalten Gases wird fiir den Fall erliutert, dass die Lewis-Zahl des kondensicrenden 
Anteils eins ist. Diese Methode kann dazu dienen vorauszusagen, ob in einem sich mischenden 
System eine Nebelbildung méglich ist. Ein experimentelles System wird beschrieben, das darauf 
angelegt ist, einige der theoretischen Schlussfolgerungen nachzupriifen. Der Apparat gestattet, 
Luft-Dampf-Gemische bekannter und regelbarer Temperatur und Konzentration herzustellen. 
Es wird tiber experimentelle Ergebnisse fiir einen turbulenten Gas-Glycerindampf-Strahl berichtet. 
Die Ergebnisse zeigen, dass betriichtlich héhere Dampfkonzentrationen als das theoretische 
Minimum notwendig sind, um eine merkliche Kondensation in einem Diisenmischer zu erhalten. 


INTRODUCTION 
CONDENSATION from a mixture of a vapour and 
a non-condensable gas is most often produced by : 
(1) cooling at the boundaries when the gas flows 
over a cold surface ; (2) adiabatic or nearly adia- 
batic expansion; or (3) mixing with another, cooler, 
gas. The first mechanism, surface cooling, 
produces a gross transfer of vapour to the wall ; 
this may be accompanied by fog formation in 


the gas if local conditions are such that supersa- 
turation is attained. This process has been 
studied by Jounstone et al. [1] who passed 
mixtures of nitrogen with vapours of sulphur, 
n-butyl alcohol and water through a_ short 
condenser section. Using an approximate theore- 
tical analysis they were able to predict in most 
instances whether or not fog would form in the 
condenser section. Their results must be regarded 
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as only semi-quantitative since no effort was 
made to measure or control the fluid mechanical 
situation. 

In the second process, adiabatic expansion, 
condensation may result when the slope of the 
isentrope on the vapour pressure vs, temperature 
plot is less than the slope of the saturation curve, 
It is the mechanism on which the Wilson cloud 
chamber is based and is of controlling importance in 
the growth of water droplets in the atmosphere [2]. 
It has been studied in some detail in high-speed 
flow processes by passing mixtures of water vapour 
and air through a converging diverging nozzle 
(see for example the review by Srever [3)}). 
Supersonic velocities are attained in the diverging 
section where both the pressure and temperature 
drop rapidly in such fashion as to produce a 
supersaturated vapour. This process is the flow 
analog of the cloud chamber. 

Condensation produced by expansion may take 
place either on foreign nuclei or as a result of self- 
nucleation (see for example Frenket [4]) depend- 
ing on the speed of the expansion process. In 
high-speed expansions in the cloud chamber or 
in nozzle flows, high values of the supersaturation 
are attained and self-nucleation dominates. In 
atmospheric expansion processes it is likely that 
foreign nuclei are most important. 

The third mechanism, mixing, is the subject 
of this paper. When a hot condensable vapour is 
mixed with a cool gas, both the vapour concentra- 
tion and temperature fall. The criterion for 
condensation is that at some point during the 
mixing process, the temperature drops sufliciently 
under the dew point to initiate stable droplet 
formation. The occurrence of 
depends on the relative rates of heat and mass 
transfer in the mixing system. It is quite possible 
that a hot vapour which would condense were 
it at the temperature of the ambient atmosphere 
will not condense during mixing. As in the case 
of expansion condensation may take place either 
on foreign nuclei or by self-nucleation depending 
on the speed on the mixing process. If mixing 
is slow condensation will take place on dust 
particles and other suspended matter in the gas. 
On the other hand, rapid mixing will produce 
high values of the supersaturation and_ self- 


condensation 


K. 


nucleation will dominate. If mixing is sufliciently 
rapid it is also possible for the system to pass 
through the two-phase region without condensa- 
tion. 

Despite the ubiquity of vapour mixing as 
a method of forming natural and artificial 
aerosols few theoretical or control experimental 
studies of the process seem to have been made. 
Among the most interesting is that of Ame.in 
and Betakov [5] who measured the size distribu- 
tion of drops formed when glycerine vapour 
was jetted into cool air. They found that the 
mean size of the drops decreased as gas velocity 
increased concluded that self nucleation 
rather than condensation on foreign nuclei was 
predominant. 

In the present paper the authors show that 
under certain conditions it is possible to trace 


and 


the path followed by the properties (composition 
and temperature) of a mixing system. In addition, 
qualitative statements can be made as to the 
effect of operating variables on particle size. 
An experimental apparatus designed to study the 
phenomenon is described and data are reported 
for the turbulent jet mixing of glycerine air 
mixtures with cool air. 


Tur Concentration Temprrarore 

A system typical of the type to be considered 
is a steam jet exhausting to the atmosphere. 
The onset of condensation depends on the con- 
centration and temperature distributions in’ the 
mixing vapours. It is possible to set up the 
equations of energy and diffusion in a number of 
ways for such a system [6]. To relate the con- 
centration and energy parameters of the mixing 
process most easily to the “ natural” thermo- 
dynamic variables of condensation a molar basis 
is chosen for the equations. In the absence of 
condensation and of thermal and pressure 
diffusion the equation of conservation of species i 
can be written in the form : 


> 


Vy Va D Vy (1) 
where the velocity v is the mean molal velocity. 
For a binary, D; can be identified with the 
binary diffusion coeflicient as defined in the 
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kinetic theory of gases. (See for example 
CharmMan and Cow1ine [7]). 
In the absence of condensation an energy 
balance on the mixture can be written in the form 
n oT VT 


(2) 


The most important assumptions involved in 
the derivation are (1) the pressure is constant 
throughout the system; (2) viscous dissipation 
can be neglected; (3) the Dufour effect (the 
interaction between diffusion and conduction) 
can be neglected ; (4) the system is a_ perfect 
gas mixture; (5) the molal heat capacities of 
all species are equal, 

Suppose a stream (or streams) of condensable 
vapour at a concentration yy and temperature 
T, mixes with gas of concentration y and 
temperature New concentration and tem- 
perature variables can be defined as: 


y — 9. 


Equations (1) and (2) become 


nw Vy’ 


n oT VT 


For the vapour jet, the boundary conditions in 
terms of the new variables are : 


y — T’ — 1 over the jet orifice 
(6) 


For steam air mixtures the Lewis number, 
DC,, is not far from one and practically 
constant over a wide range of temperatures. 
For most other gas mixtures it is of the order 
of unity. When 1 nD—kC, and 
equations (4) aud (5) are of the same form ; 
moreover the boundary conditions are identical. 
In principle with sufficiently accurate instru- 
ments it would be possible to measure v and thus 
determine the velocity field. With the same 
velocity field and identical boundary conditions 
the solutions to equations (4) and (5) must be 
the same, i.e, 


— 0 at infinity 


N LE 


This equality is independent of the geometry 

of the velocity field or of whether it is turbulent 
or laminar. 


(7) 


Tue Equation or 
Equation (7) relates the concentration and 
temperature at all points in the absence of 
condensation, The equilibrium properties of 
the system are described by an equation of the 
Clapeyron type: For the equilibrium of a pure 
liquid with its vapour in a perfect gas mixture 
this takes the form: (see for example Denricn 
diny A 
(sr 
Integrating, assuming that A is constant over 
the temperature range of interest, there results 


1 
— —{—_— — 9 
where 7, 
pressure on the system. 


boiling temperature at the total 


Tracing THE Course or CONDENSATION 

A plot of equation (9) — the saturation curve — 
is shown in Fig. 1 for glycerine vapour. The 
dilution of a vapour of initial properties y. Ty 


mole Yo 


Y= glycerine vopour conc, 


EE EE 


Fic. 1. The course of condensation in a mixing system. 
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by a vapour at yy, T is represented by the 
straight line determined by equation (7) if 
Nix = 1. Each element of vapour begins at 
concentration yy and temperature 7, and passes 
through the same series of states to y., T.. 
The time spent in each state differs depending on 
the fluid mechanics of the system. Of course in 
a turbulent flow the concentration and tempera- 
ture may momentarily If the line 
representing the locus of states does not intersect 
the saturation curve, a vapour of the given initial 


increase. 


concentration will not condense (curve A), 

For a given initial temperature and ambient 
conditions there is a minimum initial concentra- 
tion below which occur. 
This concentration is determined by the tangent 
to the saturation curve (curve B) and can be 
calculated in the following way. For simplicity 
it is assumed that yo = 0. At the point of 
tangency, y dy dT = dy, dT and 


ae, 


Also 


condensation cannot 


A;,1 1 


Solving simultancously 


RT,? — A +AT_=0 (12) 
and 
A (] —4RT_) 
4 3 
2R / A | (18) 


Only the lower root is of interest. The criterion 
for condensation becomes : 


T—T 1 
where 


A (l1—4RT_) 
2K J A | 


If the initial conditions are such that the system 
enters the saturation region nucleation and con- 
densation may take place depending on the time 
spent and supersaturation attained. The maximum 
supersaturation occurs when the system passes 
through the saturation region without condensa- 
tion. The path followed is still a straight line 


(15) 
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(curve C). The maximum value of the super- 
saturation can be caleulated by setting its 
derivative with respect to JT equal to zero, 
For y, — 0 this is given by: 


The temperature at which S,, occurs is the same 
as that at the point of tangency (equation 15). 

The initial behaviour and particle size of the 
condensing system will depend on the rate of 
mixing. If mixing is carried out slowly, con- 
densation will take place on the foreign nuclei 
which happen to be present and high values of 
the supersaturation will reached 
(curve E). Particle size will be large and concen- 
tration low. Lf mixing is rapid, the time available 
for the initial nucleation and condensation will 
be short, and it will be possible for the system 
to approach the maximum value of the supersa- 
turation where the formation of very small 
particles proceeds at a high rate by self-nucleation 
(curve D). Finally, if S,, is not large and if the 
system passes through the saturation region very 
quickly, it is possible that little condensation 
will take place at all. In any case, as mixing 
continues, — the formed must always 
become unstable and tend to evaporate since the 
ambient atmosphere is unsaturated. However, 
a meta-stable state may exist if the vapour 
pressure of the particles is low at ambient tem- 
perature. 

The effect of Np» 4 1 can be estimated quali- 
tatively. For dilute mixtures of a condensable 
vapour with air, the thermal diffusivity is usually 
greater than the molecular diffusivity of the 
vapour and N,;, > 1. Heat diffuses more 
rapidly than mass and the concentration of each 
element of gas tends to be somewhat higher at 
any temperature than the value it would have 
assumed for N;,» = 1. Henee, a gas element 
follows a path concave downwards on the y vs. T 
plot. In the absence of condensation, the path 
begins at the conditions of the nozzle and ends 
at ambient conditions. However, the series of 
states traced by each element of gas differs from 
that of its neighbour unlike the case for Np» = 1. 


never be 


aerosol 
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Fic, 2. Experimental equipment. 


EQuitpMeNT 


To study condensation by mixing, an apparatus has 
been set up which generates vapour continuously under 
controlled conditions. A diagram of the equipment is 
shown in Fig. 2. Room air is filtered and compressed in a 
j cu. ft min air pump (1). It is refiltered (2) and passed 
through a surge tank (3) and control valves (4) and (5). 
The air is dried in a Drierite column (6) which must be 
regenerated periodically by the heating system (7). 
The air is again filtered at (8). A bleed valve at (9) may 
be used to adjust to various operating conditions. Valves 
(10) and (11) and rotameters (12) and (13) are used to 
regulate and measure the flow of air to the ionizer (14) 
and vaporizer (15) respectively. The vaporizer contains 
the aerosol material. It is heated to a temperature below 
the normal boiling point by a heater (16) which is controlled 
by the system (17). The ionizer is a chamber containing 
sealed platinum electrodes. A spark is provided from power 
supplied by the system (18) which consists of a 115 V a.c. 
supply, powerstat and 10,000 V ignition transformer (19). 
Ions generated to serve as condensation nuclei can be 
combined in the vaporizer with the hot air and vapour. 

The stream from the vaporizer passes into a heating 
and mixing section consisting of a 14in. length of 1 in. 
stainless steel pipe (21) packed with small pieces of 
stainless steel screen (22). The approach section to the 
packed tube is heated by a } in. heating tape (20) controlled 
by a powerstat (27). The tube passes through a combustion 
furnace (23) and is connected by an elbow to a 25 in. 
length of 1 in. stainless steel pipe (24) which serves as the 
approach section to the nozzle (25). The nozzle is a 
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stainless steel converging section with an orifice diameter 
of } in. The approach section and nozzle are wrapped with 
a 4 in. heating tape (26) surrounded by insulating material 
(30). An iron-constantan thermocouple mounted in a 
traversing probe (31) is located 10 cm upstream from the 
mouth of the nozzle. It is used to measure the tempera- 
ture distribution in the approach section to the nozzle. 
The thermocouple leads are connected to a potentiometer 
(32) with an ice-bath junction (33) to establish the reference 
potential. 


The nozzle enters the condensation chamber (34) 
through a transite floor. The transite roof of the chamber 
is fitted with an exhaust duct which leads to a blower. 
The side walls of the chamber are perforated to allow 
passage of air, but minimize drafts. The front of the 
chamber consists of a sheet of plexiglass for visual and 
photographic observation. The condensing jet may be 
iliuminated from the side by a vertically collimated beam 
of light. 


A sampling probe (35) is located 50 diameters down- 
stream from the nozzle along the axis of the jet. During 
a run from the chamber is continuously drawn through 
the probe and then through a Sinclair Phoenix forward- 
seattering smoke photometer (36) by a blower (38). The 
photocell current is read on a logarithmic amplifier (37). 


Glycerine was chosen as the aerosol material since its 
properties are well known [10], it is easy to handle and it 
has been used in related studies. It is a “ high boiler” 
and the aerosol drops can be recovered and kept for 
examination. The glycerine employed was reagent grade. 
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Prockpure 

A known weight of glycerine was placed in 
the vaporizer (15). With inlet valves (10) and 
(11) shut the vaporizer was heated to approxi- 
mately the desired operating temperature. 
During this period the entire system was heated 
to prevent condensation in the line during the 
start of a run. The loss of glycerine in the boiler 
during this time was negligible. Next, either 
or both of the valves (10) and (11) were opened 
and the flow rates were maintained constant, 
The temperature in the vaporizer reached a 
constant value within 3-4 min. The generation 
system was allowed to operate in this fashion 
for from 4 7 hr. The average concentration of 
vapour in the efflux was determined from the 
weight loss of glycerine and the known air-flow 


rate. 

The presence of fog in the jet was detected by 
the photometer (36). Before beginning a run 
the chamber was purged with room air and the 
meter (37) adjusted to zero. This adjustment 
showed little fluctuation from day to day. The 
vapour and air from the boiler were heated in 


the packed tube and entrance section until no 
condensation could be detected. A rough deter- 
mination of the initial condensation point was 
made by rapidly lowering the nozzle temperature, 
T,, until condensation was noted. The gas stream 
was then reheated until the meter reading re- 
turned to zero and a careful determination made 
by lowering the temperature slowly over a period 
of 30 min or more. At some point the meter 
reading rose sharply to values of 2-5 correspond- 
ing approximately to a 10®- 10° increase in mass 
concentration of condensed matter [9]. This 
point was defined as the onset of condensation, 
A check was made by slowly raising the tempera- 
ture until the meter reading returned to zero, 
Within experimental error, the same point was 
found. At the condensation point, it was possible 
to observe visually occasional puffs of fog in the 


chamber. 


Discussion 
When a vapour stream is mixed with cool air 
there are ranges of initial concentration and 
temperature which do not lead to condensation. 


This was observed in a preliminary investigation 
in which photographs were made of the condens- 
ing jet for various initial concentrations of 
glycerine in the range yy — 0-003 per cent to 
4% per cent. The initial temperature, 
Ty. was 190°C and the gas velocity was approxi- 
mately 1000 em sec. Figs. 3 (a) (e) show the 
effect of decreasing concentration of the position 
of the condensation zone. As y, decreased, this 
zone tended to blow off the nozzle until it dis- 
appeared completely. 


00028 00030 


Fic. 4. Limiting conditions at nozzle for condensation 
in a turbulent jet: nozzle velocity ~ 1000 em sec. 
A with ions. © without ions, 


The experimental programme was concerned 
largely with investigating the limits of vapour 
concentration and temperature at the nozzle 
which lead to detectable condensation. The limits 
were measured for glycerine air streams both with 
and without the introduction of ions. Glycerine 
concentrations at the nozzle ranged from 0-038 
to 0-273 mole per cent and temperatures from 
197°C to 292°C. In this study it was not intended 
to vary the gas velocity at the exit of the nozzle 
but values did range from 900 to 1380 em sec. 
The results are presented in Fig. 4. The region 
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above the lines represents vapour temperatures 
and concentrations at the for which 
condensation was found in the jet. No detectable 
condensation took place for nozzle conditions 
corresponding to points in the region below the 
lines. The mass concentration of the fog increased 
sharply (about three or four orders of magnitude) 
when the ions were introduced, No measurement 
of ion concentration was made but the air velocity 
through the ionizer and the applied potential 
were kept constant in all runs. 

Also shown in Fig. 4 are the saturation curve 
for glycerine and the tangent curve (line B in 
Fig. 1). The saturation curve is based on the 
data of SrrpMaN [11] assuming a total pressure 
of Latm. An independent determination of the 
vapour pressure of glycerine has been made by 
Fitosoro et al. [12]. Their results agree with 
those of SrepMan. The results of Zit’ BeERMAN- 
GRANOVSKAYA as reported in [10], (p. 266) do not 
agree with SrepmMan’s measurements. The 
tangent curve is drawn for — 28°C, 0, 
and a total pressure of 1 atm. 

If the Lewis number were unity for glycerine - 
air mixtures, the tangent line would represent 


nozzle 


the lowest position along which the experimental 
curves could fall, and Fig. 4 would indicate that 
concentrations considerably in excess of the 
theoretical minimum values are necessary for 
Actually, for glycerine - 
air mixtures NV,» ~ 3 based on a diffusion coetli- 
cient of 0-0695 em see (20°C, Latm) calculated 


detectable condensation. 


from the Arnold equation. The path followed 
by a gas element on the y vs. T plot is concave 
downwards since, as discussed above, cooling 
occurs more rapidly than mixing. Hence even 
lower concentrations should be necessary at 


the nozzle to produce condensation for N;» = 3 
than for N,;, — 1. This strengthens the con- 
clusion that very high supersaturations must 
be attained in a rapid mixing process for detect- 
able condensation to occur. In the presence of a 
large amount of inert (air) the growth of nuclei, 
whether foreign or self-generated, probably de- 
pends on a process of local diffusion. This 
process is too slow to permit the growth of nuclei 
to a size large enough for detection by the photo- 
meter unless the local supersaturation is very 
high. In the absence of condensation, it is pro- 
bable that most of the glycerine is dissipated 
molecularly or as small clusters which soon re- 
evaporate as mixing proceeds. 
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NOTATION 
molal heat capacity 


diffusion coeflicient 
thermal conductivity 
molal density 
Lewis number = k/n DC, 
gas constant 
supersaturation 
- time 
temperature 
molal gas velocity 
concentration of condensible vapour, mole 
fraction 
- molal latent heat of vaporization 


inlet conditions 
ambient conditions 
tangency conditions 
- maximum value 
- boiling 
8 — saturation curve 
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On laminar boundary layer flows with a rapid homogeneous 
chemical reaction 
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Abstract—A theoretical analysis is presented of a mass transfer process which is coupled with 
an instantaneous homogeneous chemical reaction between two substances. It is shown how 
the appropriate diffusion and energy equations for laminar boundary layer flows past arbitrary 
surfaces may be solved in closed form by making use of an asymptotic technique. Of particular 
interest is the derived expression for the rate of mass transfer, of one of the reacting species from 
the surface, which has the simple form 


Nu = Nuy C (x, A) 
Nuyg is the local Nusselt number in the absence of a chemical reaction, and C (a, A) is a correction 


factor, independent of the surface geometry, which can be calculated exactly as a function only 
of two fluid parameters « and A. 


Résumé Présentation de lanalyse théorique d'un processus de transfert de matriére couplé 
avee une réaction chimique homogéne instantanée entre deux substances. On montre comment 
résoudre par une méthode de tangentes une expression implicite des équations convenables pour la 
diffusion et Pénergie pour écoulement Laminaire de la couche limite au contact de surfaces arbitraires. 
Pour la vitesse de transfert d'un des réactifs au travers de la surface, on peut dériver une expression 
particuliérement intéressante de forme simplitiée 
Nu Nuy C (x, A) 

Nu, étant le nombre de Nusselt local dans Pabsence de réaction chimique et C (x, A) est un facteur 
de correction, indépendant de la géomeétrie de la surface, pouvant étre calculée calculée en fonction 
de deux paramétres du fluide « et A. 


Zusammenfassung -Eine theoretische Analyse cines Stoffiibergans wird dargelegt, der mit einer 
rasch verlaufenden chemischen Reaktion zwischen zwei Substanzen gekoppelt ist. Es wird 
gezeigt, wie geeignete Diffusions- und Energiegleichungen fiir laminar, tiber beliebige Oberflachen 
fliessende Grenzschichten dadurch in geschlossener Form gelést werden kénnen, dass man eine 
Asymptotentechnik benutzt. Von besonderem Interesse ist der Ausdruck, der fiir die Massenii- 
bergangsgeschwindigkeit ciner der beiden Reaktionspartner von der Obertliche weg abgeleitet ist. 
Er hat die folgende Form 
Nu = Nuy C (a, A). 

Nuy ist die Ortliche Nusselt-Zahl unter Ausschluss der chemischen Reaktion. C(x, A) ist ein 
Korrekturfaktor, der von der Gestalt der Obertliche unabhianging ist und der als Funktion der 
beiden Parameter des fliessenden Mediums « und A genau berechnet werden kann. 


INTRODUCTION fluid stream will be affected by the presence of a 
Fivrm mechanical phenomena involving the chemical reaction, either on the surface itself 
simultaneous transfer of matter with a chemical — or in the bulk of the fluid. In the present theore- 
reaction have received a considerable amount of — tical analysis we shall deal with a special problem 
attention in the last few years [1, 2, 3, 4]. One in this broad field. We shall consider the effect 
of the reasons for this is that it is of some impor- on the rate of mass transfer of a very rapid 
tance in a number of chemical engineering opera- homogeneous reaction between two chemical 
tions to be able to determine to what extent the — species, and show how the solution of FrrEDLANDER 
rate of mass transfer from a surface to a moving and Lirr [5] to the equations of diffusion and 
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energy, for such a reaction occurring in the 
laminar boundary layer on a flat plate, may be 
generalized to systems with an arbitrary surface 
geometry. 

Let us visualize a typical physical set-up for 
such a process. A fluid containing a reactant B 
flows past the arbitrary surface, shown schemati- 
cally in Fig. 1, which generates the other reacting 
species A. As an example, A might be an organic 
acid which is slowly dissolving from the solid 
surface and instantaneously reacts with a liquid 
stream containing a base, Or, alternatively, if 
the convective fluid is a gas, A might be a vapour 
generated from the surface either by evaporation 
or sublimation. In any case, since the reaction 
between 4 and B is assumed instantaneous, it is 
clear that the flow field will be divided into two 
distinet regions which, as is illustrated in Fig. 1, 
are separated by a reaction zone of infinitesimal 


Fic. 1. A schematic diagram of the physical set-up. 


In order now to achieve our main objective, 
which is to ascertain the effect of the chemical 
reaction on the rate of mass transfer of substance 
A, we shall limit ourselves to a hydrodynamic 
system, the laminar boundary layer flow, which 
has the unique property of being well enough 
detined to admit of a theoretical treatment, 
while at the same time being realized often 
enough in practice to be of more than academic 
interest 


Basic Equations axnp Turik 


The process under consideration is governed 
by the diffusion equations for convective mass 
transfer of the two species A and B which, for 


laminar boundary layer flows, become [5, 6] 


with the boundary conditions : 


= €,(0) 
= 0 

=0 

Ca = Cy (®) 
Cy(®) 
= 0 


aty, 0 

ata, = 0 

aty, = 

atry, 0 


(3) 


where ¥p (7,) denotes the position of the reaction 
zone. In addition, if the chemical reaction takes 
place according to the stoichiometric equation 


aA ~ bB -» products (4) 


it is clear that a material balance around the re- 
action zone will vield the last boundary condition, 


a hb 
Seg Sen 


In what follows, we shall suppose that the two 
velocity components u, and v,, to be substituted 
in equations (1) and (2), may be obtained from 
the solution of the well-known laminar boundary 
layer equations for the transfer of momentum 
|6, 7). It is realized of course that it is not quite 
correct, strictly speaking, to assume that the 
hydrodynamic equations may be solved indepen- 
dently from the mass transfer equations, for not 
only are the fluid properties in general dependent 


(5) 


on composition but, as a direct consequence of the 
transfer of material at the surface, the normal 
component of the interfacial velocity is no longer 
zero, It is known, however, that the effects just 
mentioned are usually of secondary importance 
only and will therefore be disregarded for the 
purposes of this study. 

Now, Frieotanper and Lrrr [5] have already 
shown that equations (1) and (2) may be solved 
analytically, by a similarity transformation, only 
for flows past a wedge-type surface, of which 
the flat plate is but a special case. It is the purpose 
of the present analysis to remove this restriction 
and to obtain a solution to such a physical system 


1 d% 
+0, — (1) 
wy 
and 
d 
13 
26 
> i 
— 
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which will apply to two-dimensional surfaces of 
arbitrary geometry. This will be accomplished 
in the following manner: It was noted some 
time ago by Face and Fatkner [8] and proved 
rigorously by Morcan and Warner [9] for the 
analogous heat transfer problem that when the 
Schmidt number Se is large it is permissible to 
replace the exact velocity components u, and x, 
in the diffusion equation by their asymptotic 
form near the surface. Furthermore, it was 
established by Licuruins. [10] in an often quoted 
paper that this asymptotic development, which 
strictly speaking becomes exact only when 
Se is usually of acceptable accuracy even 
if Se is around unity. This means therefore that, 
since for most gases Se ~ 1 and for all liquids 
Se » 1, we may confine our attention to the 
asymptotic solution of equations (1) and (2) 
when both Se, © and Sey ©. For, in view 
of Licurnit’s results, we should expect that any 
conclusions derived from such an asymptotic 
analysis of our problem would also apply to 
essentially all possible cases of practical interest. 

Thus, when both Se, » 1 and Sep, » 1, the 
physical process of mass transfer and chemical 
reaction will take place in a region very close to 
the surface in which u, and vy, may be represented 
respectively by 


m= Blam 
and 
2 (6) 
y*, dp 
vy — — 
2 dr, 
The function 8(7,) is proportional to the shear 
stress at the surface which is either known or 
may be calculated for a large number of geometries 
by standard methods [7, 11, 12]. It is possible 
to show now that equations (1), (2), (3). (5) and 
(6) can be solved analytically by a similarity 
transformation. Let 


(n) 


=O, (»); Cy 


(0) 
where 
7 =m (Se,)'* 1/3 (7) 
| VB de 


which reduces our system of equations to 


+ =0 
and 


with the boundary conditions : 


1, O4(nR) = 9 
On (np) 0, (@) 1 (9) 


and 
AOn(n,) 


nx denotes the position of the reaction zone, 
defined by equation (7) with y, = yp (#,), and « 
and A are two known fluid parameters given by 

Sc 

and A= omy Cn 

b Sep (0) 
Equations (8) and (9) may now be solved by ele- 
mentary means. Thus 


: Fexp(— 2)'dz 
. 


| 2 exp (— 2) dz 
. 


x 


| >* exp (— 2) dz 

where »,. Which is a function only of A and «, 
must be obtained from 


exp (— 2) dz 


— exp [(2 — 1) nal (12) 


| = exp (— 2z) dz 


Since the integrals appearing in equations (10), 
(11) and (12) are expressible in terms of the in- 
complete gamma-function, already tabulated by 
Pearson [13], the numerical evaluation of 
0, and @, can be performed quite readily. Of 
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0,=1- for 0 <» (10) 

| exp(—z)dz 
and 
On for @ (11) 
| 
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particular interest is the rate is mass transfer 
of A. If 


then, in view of equations (7) and (10) 


rs) 
Nuy VR 
exp (— 2) dz 


where Nu, is the Nusselt number in the absence 
of any chemical reaction (A — 0 or yn, » %) and 
is equal to 


vp va (14) 


3) (| vex) 


which is incidentally identical to the well-known 
Licuruitt formula [10] for the analogous heat 
transfer problem. 

The principal result of our analysis may now 
easily be deduced from equation (13), for it is 
clear that, since yn, according to equation (12) is 
determined solely by the numerical values of « 
and A, the effect of the chemical reaction on the 
rate of mass transfer of A can be represented 
by a simple expression of the form 


Nu = C Nu, (15) 


where the correction factor C which may be 
evaluated exactly from equations (12) and (13) 
is of a function only of the two dimensionless 
groups x and A and is completely independent of 
the surface geometry. We have succeeded in 
establishing, therefore, that equation (15) which 
had already obtained by Frieptanper and Lrrr 
[5] for the flat plate only, is valid for any two- 
dimensional surface. 

Before concluding this section we shall briefly 
derive some useful properties of the function 
C (a, A). 


(1) A closed-form expression for C (x, A) may 
easily be obtained when x — 1, that is, 
where Se, = Sey. Thus it is obvious from 
equations (12) and (13) that under this 
condition 


C(x Ay=1+A 


which, as can be verified from the basic 
differential equations themselves, is valid 
even where Se, < 1. 


When A > © one would certainly expect 
that the rate of mass transfer of A from 
the surface would be completely controlled 
by the rate of transfer of B from the 
bulk to the reaction zone. This is borne 
out by our results. Thus, as A -+ », 
0 and therefore 


C (a, A) >» Ax'® as A + ow (16) 


The general behaviour of the function 
C (x, A) between its two asymptotes, 1 and 
Ax! *, is depicted in Fig. 2 for x = 10, 1 
and 0-1, while »,, the predicted position 
of the reaction zone, is shown in Fig. 3. 


Fig. 2. The general behaviour of the function C (a, A). 


20 


02 04 06 
Ae 
Fig. 3. The position of the reaction zone as a function 
of the two fluid parameters a and A. 


Tur or Tue Enercy Equation 


Our general analysis can finally be extended 
to non-isothermal processes and may be used to 
predict the temperature rise in the reaction zone. 
This information is usually of some importance 
whenever chemical reactions take place and where 
large heat effects must be avoided. 


VOL. 
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For two-dimensional laminar boundary layer 
flows the temperature distribution is governed 
by the energy equation [6, 7] 


1 YT 


Pr dy,’ 


17) 


If again we make use of our asymptotic develop- 
ment with Se, » 1 and Se, » 1, we may replace 
u, and v, by equation (6) since the chemical 
reaction will occur in a region very close to the 
surface. Furthermore, if we let 


T = T(») 


where, as before, » is given by equation (7), we 
can shown that 


dT 
dr? dy (18) 


We shall now solve equation (18) subject to 
the boundary conditions : 


(a) T=T, at 0 and » = x 
AH, Pr do 

” Se, (i, 


dT 


which is obtained from a heat balance around the 


» 
reaction zone. Thus, if A, = (0), we 
pC, Se4 
can readily verify that 
(a) For » <p 
T—T, _ exp(% — 
[= exp (— 2) dz a" 
(20) 


exp (— 2) dz° 
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r = ~ 1) or 


exp (— 2) dz 
j fe exp(—2)dz (20a) 
2/8 exp (— 2) dz 


0 


Of particular interest is the maximum value of 
the temperature, 7. which is obtained from 
equation (20) by letting » = nz. Under special 
conditions two simple expressions can be derived : 


(a) When x, = 1, which is usually obeyed 
approximately by most gases, 


2 = (21) 


where C is given by equations (13) and (15). 


(b) When A -+ «, that is, when the process is 
controlled by the rate of transfer of B from the 
bulk to the reaction zone, 


(22) 


It is, moreover, easy to show, from an analysis 
of the basic differential equations, that this is a 
more general result which is not restricted to 
the simple system considered in this work. 


CONCLUSIONS 

In this paper we have presented a theoretical 
analysis of a mass transfer process which is 
coupled with an instantaneous homogeneous 
chemical reaction between two substances. We 
have shown how it is possible to obtain a closed- 
form solution to the appropriate equations of 
diffusion and energy in laminar boundary layer 
flows past two-dimensional surfaces of arbitrary 
geometry, by making use of an asymptotic 
development which, although strictly speaking 
exact only when the Schmidt numbers of both 
diffusing species approach infinity, is of acceptable 
accuracy even when Sc ~ 1. As a consequence 
the results of our present work should apply 
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to most physical systems of interest under laminar 
boundary layer flow conditions. 

The principal conclusion of this analysis 
concerns the way in which the rate of mass 
transfer of one of the reacting species from the 
surface is influenced by the presence of the 
chemical reaction. We have been able to show 
that, if Nu is the local Nusselt number and Nuy the 
the Nusselt number in the absence of a chemical 
reaction, then 
Nu — Nuy (x, A) 


where the correction factor C (x, A), which may be 
evaluated exactly from equations (12) and (13), is 
a function only of the two fluid parameters « 
and A, and is independent of the surface 
geometry. 
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a,b — stoichiometric coefficients of the chemical 
reaction 
A, B = the two reacting components 4 and B 
Cy = specific heat 
C4 = concentration of component 1 
¢, (0) = concentration of A in the fluid next to the 
surface 
concentration of component B 
€y (©) = concentration of B in the bulk 
C (a, A) = a function defined by equations (13) and (15) 
D — diffusion coeflicient 
k = thermal conductivity 
L = characteristic length 
T = temperature of the fluid 
T,. = temperature in the bulk of the fluid 
UU, = characteristic velocity of the stream 
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a2 = distance along the surface from the leading 


= a-component of the velocity (4, = u/U,) 
v = y-component of the velocity 


Rev 


edge a, = L) 


=~ distance normal to the surface 


(yy = Vv Rey L) 
the distance of the reaction zone from the 
surface 
dimensionless fluid parameter defined by 
a = Sep/Se4 
dimensionless fluid parameter defined by 
a) Pr/Se 
1 toh 
Vv Re pl 


dimensionless concentration @ 4 ¢4/¢4 (0) 


dimensionless concentration Op 
similarity variable defined by equation (7) 
the position of the reaction zone [defined by 
equation (7) with = 
dimensionless fluid parameter ; 

ase (®) 

bSe (0) 

fluid parameter with units of temperature ; 


A 


p 4 
Viscosity 


kinematic viscosity 

mass density 

local shear stress at the surface 

the heat of reaction per unit mole 

(All, > 0 for an exothermic reaction) 

the local Nusselt for the transfer of 4 ; 

Nu L (00.4) 29 y 

the local Nusselt number in the absence of 
a reaction (A = 0) 


», the Reynolds number 


C, the Prandtl number 


v/ D, the Schmidt number 
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Abstract Numerical solutions are presented for differential equations representing the iso- 
thermal steady flow tubular reactors with varied degrees of backmixing of the reactants in the 
axial direction and with reactions with order other than zero and first. 


Résumé—Les auteurs présentent des solutions numériques pour les équations différentielles 
représentant écoulement stationnaire isotherme dans des réacteurs tubulaires présentant des 
degrés divers de mélange en retour des réactifs dans le sens axial et pour des réactions dont 
Vordre est différent de zéro ou de un. 


Zusammenfassung. Es werden numerische Losungen von Differentialgleichungen 
isotherme kontinuierliche Rohrreaktoren angegeben, wobei cine verschieden starke Riickver- 
mischung der Reaktanden in axialer Richtung beriicksichtigt wird. Diese Loésungen gelten 
fiir alle Reaktionsordnungen ausser der nullten und der ersten, 


INTRODUCTION 2U — UR, =0 


A watrentat balance in a tubular flow reactor dy? dy 


operated under isothermal steady state conditions 


ul 


yields the following equation ; where U oi ead &, = — 
» 


u 

dl 
at dz r(c)=0 (1) The boundary conditions are transformed to 
When the reaction is of n-th order involving 
one reactant or several reactants with equivalent 
concentration, the rate of reaction r (e) is repre- 
sented as: 1, 


0 : — 1) 


For the zero and the first-order reactions equation 
(3) yields linear ordinary differential equations 
0. dC ic and analytical solutions are available for them 
(1, 2]. 
dc For chemical reactions other than zero and first 
=L, -—* 0 order equation (3) yields. non-linear second-order 
ordinary differential equations. Anayltical solu- 
The following transformations of the variables — tions for them are not available. Because of 
lead to the dimensionless form of differential — their importance for reactor design numerical 
equation which contains two parameters; one, U. solutions were attempted by the authors. 
being associated with the ratio of flow and diffu- 


sion velocities and the other, R,, with the ratio NumenicaL SOLUTION 
of reaction and tlow velocities. 


The boundary conditions are as follows [1]: 


The truncated finite-difference representation 
Cc of the differential equation given by equation (3) 


Cy L yields [3] 
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Y,-,(1 + Uh) +- — Uh) - 
—~UR,Y/{h? —0 (5) 


where / is the constant interval between pivotal 
points and / refers to i-th pivotal points. 

The boundary conditions given by equation (4) 
may be represented as [3] 


g=1, ¥, (6) 


7=0, Y,— + 4hU=0 (7) 


where p is the number of pivotal points. 
Several criteria were used to insure the precision 
of the numerical solution to the finite difference 
equation. Only slight improvement precision 
was obtained by increasing the number of in- 
crements used in the solution of the finite equation 
when the number of increments exceeded fifty. 
The boundary condition at » 
by equation (6) was used as the starting point in 
the solution of the finite difference equation. The 
values of Y at the pivotal points were then 
obtained by means of equation (5). The boundary 
0, represented by equation (7), 


1 as represented 


condition at » 
was used to check the calculation. 
value of Y at » — 1 was systematically changed 
by incremental values and the residue of equation 
(7) computed. The assumed value of Y leading 


An assumed 


to the smallest residue gives the closest solution 
to the differential equation. Smaller increments 
of Y at » = 1 insure residues approaching zero, 
ie. greater precision. In the data presented here 
the final increment used was 0-1 per cent of the 
value of Y at » = 1-0, and the number of pivotal 
points used set at 50. Further increase in precision 
could be obtained but would not alter the curves 
presented here. 


Resvuuts Discussion 


Results of the numerical computations des- 
cribed in the previous section are plotted in 
Figs. 1-14 for }-, }-, second- and third-order 
reactions. (Similar diagrams for any combination 
of n, R, and U may be obtained from the authors). 
Computations have been made for the case of 
second-order reactions by Levensriet and 
Biscnorr [4]. The effect of axial diffusion on the 
concentration profile and degree of conversion of 
reactant is clearly illustrated in those figures. 
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Fic. 1. Concentration profile of reactant in reactor 
for j-order reaction when KR, O2 KR, = OF 
v = Parameter. 


The figures also show that the backmixing has 
profound effect on the concentration profile in 
the reactor in all the cases. It has, however, little 
or negligible effect on over-all conversion when 
the order of reactions is small, (theoretically it 
has no effect for zero-order reaction) especially 
when is small. 
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ric. 2. Concentration profile of reactant in reactor 
for j-order reaction where R, = 1 n=} R, = 1 
v = Parameter. 
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When the axial diffusion of reactant is negligible When the axial diffusion is intinitely large and 


equation (3) reduces to: the concentration of the reactant is uniform in the 
dyn (10) 


Consequently, for » other than 1: 


100 
Y=j1-— (9) 
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090 
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. Fic. 5. Concentration profile of reactant in reactor 
for j-order reaction when R, = 1 n=} 
v = Parameter. 
Fic. 3. Over-all conversion as a fraction of reactant 
fed for }j-order reaction, n = | K, = Parameter. 
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for }-order reaction when R, = 02 n= } Fic. 6. Over-all conversion as a fraction of reactant 
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Fic. 7. Concentration profile of reactant in reactor 
n=2 


for second-order reaction when R, = 02 = 
= Parameter. 
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Fic. 9. Concentration profile of reactant in reactor 

for second-order reaction when R, = n=2 
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Fic. 8. Concentration profile of reactant in reactor 
2 Fic. 10. Over-all conversion as a fraction of reactant 


for second-order reaction when RK, = n= . 0. 
fed second-order reaction. n = 2 RK, = Parameter. 


v = Parameter. 
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Fic. 11. Concentration profile of reactant in reactor 
for this order reaction when R, = 02 n= 3 
v = Parameter. 
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Fic. 12. Concentration profile of reactant in reactor 
for this order reaction when R, 1 n= 3 
Parameter. 
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Equations (9) and (10) were used as criteria for 
validity of the numerical analysis employed. By 
increasing the values of U the concentration 
profile approaches asymptotically to that repre- 
sented by equation (9), while by reducing the 
values of U it approaches the concentration 
profile represented by equation (10). 

The results of this work show that the extra- 
polation or interpolation of concentration profile 
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Fic. 13. Concentration profile of reactant in reactor 
for third-order reaction when R, = 10 n = 2 
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for the various order of reactions with the various 
degrees of backmixing is difficult. 
separate numerical solutions as illustrated in this 
article would be necessary for each different 
condition of reactor operation. 


Therefore 
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Fig. 14. Over-all conversion as a fraction of reactant 
fed third-order reaction. n = 3 RK, = Parameter. 


Chem. Engng. Sci. Vol. 13, No. 2. October 1960. 


+— 
| 
095 + 
094 
093 100 
~ 
092 | | 
om oon 
090 ac | 
| 
| oeo} 4 4+—+— 
| 
or 
= 
MOL. 
13 \ > > 
96 \ 
050 + + + 
100 S 
4 
080 + 
075 _ 
070 
os 10 
| 
- 
Kos + t+ + — 
| | | | 
Ss / 
| | Ss 
| 
a 
U 
67 


Cc = 
— inlet concentration of reactant 
= axial diffusivity of reactant 

= mean value of E 
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NoTATION 
concentration of reactant 


constant interval 
rate constant of n-th order reaction 


= length of reactor 
= order of reaction 


number of pivotal points in finite difference 
calculation 


2k, L Cc,” 1 
u 


mean linear velocity of reacting mixture 
al 
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Zusammenfassung—Fiir die Messung der effektiven Diffusion und effektiven Wirmeleitung 
in radialer Richtung in einem festen oder fluidisierten Bett ist eine Anordnung geeignet, bei 


der der radiale Austausch zwischen zwei konzentrischen Bereichen der Stromung erfolgt. Die 
fiir die Auswertung erforderliche Integration der Differentialgleichung wird durchgefiihrt unter 
Beriicksichtigung von Wandreflexion und longitudinaler Diffusion. Die numerische Ermittlung 
der Ldsung bezieht sich auf einen Bereich, in dem die longitudinale Diffusion vernachlissigt 


werden darf. 


Abstract——Effective diffusion as well as effective thermal conductivity in transverse direction 
in a fixed or fluidized bed can be measured in a column where radial exchange takes place between 


two concentric areas of the flow field. In this article an integration of the underlying differential 
equation is given, allowing for reflexion at the tube wall and for longitudinal diffusion. The 
. solution is evaluated numerically for a range of conditions where longitudinal diffusion may be 


neglected. 


Résumé—On peut mesurer la diffusion effective et la conductivité thermique effective dans 


le sens transversal d‘un lit fixe ou turbulent, dans une colonne lorsque léchange radial s’effectue 
entre deux zones concentriques du champ d’écoulement. L‘article fournit une intégration de 
léquation différenticlle, permettant de tenir compte de la reflection sur la paroi et de la diffusion 
longitudinale. La solution numérique est évaluée numériquement pour un intervalle de conditions 


Zurn Messung der diffusionsartigen Querausbrei- 
tung (transversale effektive Diffusion) in Fillkor- 
persaulen oder verwandten Medien leitet man 
gewohnlich einen geeigneten Indikator zentrisch 
und kontinuierlich in die Strémung ein. Zur 
einfachen Auswertung wird die Quelle als punkt- 
formig angenommen. Aus der Konzentrations- 
verteilung in einem  stromabwirts 
Querschnitt gewinnt man die modellmassige 
Kenngrésse fiir die Querausbreitung. Hierbei 
kann mit Hilfe einer Lésung von KLINKENBERG, 
Krasensernk und Lauwerier [1] auch der 
Einfluss einer Langsausbreitung beriicksichtigt 
werden, 

Bei einer Untersuchung tiber Querausbreitung, 
iiber die an anderer Stelle berichtet wird [2], 
haben wir andere Randbedingungen gewahlt. 
Die indizierte Flissigkeit wird durch eine Roéhre 


yelegenen 


lorsque la diffusion longitudinale peut étre négligée. 


2A” CC, | 2A 
c=0 
: 
Ans. 1. Schematischer Lingsschnitt durch die mit 


Fillkérpern gefiillte Siule. 


(Durchmesser 2R) zentrisch zur Hauptstromung 
(Saulendurchmesser 24) cingeleitet (Abb. 1). 
Die Strémungsgeschwindigkeit w ist in beiden 
konzentrischen Bereichen iiber die gesamte 
Saulenlinge gleich. Im Abstand L von der 
Aufgabestelle wird durch ein ebenfalls zentrisches 
Rohr gleichen Durchmessers 2R die Messfliis- 
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sigkeit entnommen und aus ihrem mittleren 
Gehalt an Indikator die Kenngrésse der Queraus- 
breitung berechnet. Diese Methode vermeidet 
die mit der Realisierung einer Punktquelle 
verbundenen Schwierigkeiten; ein weiterer 
Vorteil ist die Mittelung tiber alle Azimutrich- 
tungen und die Entnahme der Messflissigkeit im 
Innern der nahezu ungestérten Packung. Die 
genannten Randbedingungen eignen sich nicht 
nur zur Messung der effektiven Diffusion in 
Festbetten fluidisierten Betten, sondern 
auch zur Messung der effektiven Warmeleit- 
fahigkeit. Wir méchten deshalb die zur Auswer- 
tung erforderliche Rechnung und ihre Ableitung 
hier gesondert mitteilen. 

Wir bezeichnen die Koeflizienten der diffusions- 
artigen Ausbreitung (Mischungskoeflizienten) mit 
M,, und M,. Da die Ausbreitung im allgemeinen 
nicht isotrop erfolgt, wird zwischen dem _ trans- 
versalen und dem longitudinalen Mischungsko- 
eflizienten unterschieden. Die Langsausbreitung 
ist jedoch in manchen Fallen wegen des geringen 
Konzentrationsgradienten in Achsenrichtung zu 
vernachlassigen. Wann dies erlaubt ist, muss 
eine Uberschlagsrechnung zeigen, die an der 
Funktion fir die  Konzentrationsverteilung 
durchzufiihren ist (Gleichung 3). 

Aus dem Elementargesetz det  anisotropen 
diffusiven Ausbreitung, j M-Ve + ev, wobei 
M ein Tensor (zweiter Stufe) mit der Matrix 


M,0 0 
» 
o M, 


ist, folgt als Kontinuitatsgleichung fir den sich 
einstellenden stationaren Zustand die Differen- 
tialgleichung 


Ye  M, we 


wv? M, dz 
(Vergl. die Erklarung der Zeichen am Schluss der 
Arbeit). 


Es wurde M, 
angenommen. 


und M, M, 


Als Randwertbedingungen werden eingefiihrt : 


dc 


r= A: 
or 


= 0 (Wandreflexion) 


r<R: we, 


dc 
=we— M, 
dz 


R<r<aA: = w — M,% 
dz 


d.h. die Einlassebene wird als undurchlassig 
gegentiber diffusiver Rickmischung angenommen. 
In der Messebene (z — L) sei der mit der z-Rich- 
tung abnehmende axiale Konzentrationsgradient 
schon so klein, dass ein Einfluss des Saulenauslaufs 
infolge diffusiver Rickmischung zur Messebene 
vernachlassigt werden kann; die Saule wird also 
in Strémungsrichtung als unendlich angenom- 
men, 
Nach Einfihrung dimensionsloser Gréssen* 


dw dw 
in denen d eine fiir die Packung charakteristische 


Lange. z.B. Teilchengrésse, ist, lauten die Differen- 
tialgleichung und ihre Randbedingungen : 


S(p) =y- 


R 


Bodt 
Ein Ansatz der Form y — P(p)- fihrt auf 


die vertragliche, tberall endliche Partikularlé- 
sung: 


do (p A,) 


exp (Bo, |(1 + - 
J,(p A,) — Besselfunktion der Ordnung Null. 


*Auf Vorschlag von D. W. van Krevecen [3] wird die 
Kennzahl der diffusiven Ausbreitung nach Bodenstein 
benannt; sie entspricht der fiir Temperaturleitung 
definierten Pécletzahl. 


y=C, 


VOL, 
13 
Bo, My _ Bo, —0 (1) 
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Die Lésung des Randwertproblems werde als 
eine Summation tiber solche Partikularlésungen 
angesetzt (Fouriersche Methode) : 


ZC, Jo (p A,) 


Aus der Randbedingung 


k* Bo,, Bo, 


A? 


1/2 
folgt J, (A,) = 0. 


Die Summation ist also tiber alle Wurzeln der 
Gleichung J, (A) = 0 zur erstrecken : 


No 3.8317 ; Ay = 7,0156 ; 


A; = 


= 0; Ay = 
10,1735 u.s.f. 


Zur Bestimmung der C, benutzt man die Rand- 

bedingung der Gleichungen (2b). Die Funktionen 

u, = Jy (pA,) mit J, = 0 bilden einen 

vollstandigen Satz von Orthogonalfunktionen. 

Die C, ergeben sich mit Hilfe der Orthogonalei- 


genschaft | pu,u, dp = 0,(v Es sei weiter 


0 
ein beliebiger Wert unter den natiirlichen Zahlen 
v. Dann errechnet sich schrittweise (vgl. Defini- 
tion von f(p) in Gleichung (2b)) : 


S(p) p Jo A,) = ZC, r) Jo(pA,) 


pS (p) Jo (p A,) dp = 


c|(t+ 8] fon (p A,) dp. 


R/A 


+ Bo, Bo,)]'* + 4 
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Die Konzentrationsverteilung folgt somit der 
Funktion 

R J, (A, RA) 
exp | — Bo,[[} + (* Bo, 

Fiir den Term v = 0 ist hier und in der Folge zu 
beachten, dass 


2 


y= (A, p) 


(3) 


J, (0) = lund ilt. 


= 


Durch jeden Querschnitt der Saule muss wegen der 
Stationaritaét des Vorganges die gleiche Indikator- 
menge fliessen. Dass diese Bedingung erfiillt ist, 
zeigt eine Integration itiber den Flussvektor 
jz = we — M,de/22: 


[iar = A? Bo p dp. 
(F) 
Es ist 
1 
| w dp 
— CBo,[[t + (A? Bo, — 
[4 + (A,2/k? Bo,, + § 
1 
x | Jo (A, p) p dp, 


R 
A 
v=0 


J, (A,R A) 
A, (A,) 


1 
_ ky 
A 
v=0 


J, (A, R/A) 
A, 


Bo, [[} + (A,2/k* Bo,, Bo,)]'’? — 
[4 + (A,2/k* Bo,, + 


alt 


1 oY 
d¢ 


A? 


1 
Das Integral ) dp nimmt far A,#A,=0 


den Wert Null an wegen 
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1 
J, (A) 0A, 0. 


1 
| (A, p) p dp = 
0 


1 
Im Fall A, = Ay = 0 folgt eae = 


Die Integration tiber den konvektiven Transport 
ergibt also fiir sich allein einen von der Achsen- 
koordinate z unabhaingigen Wert, da fir Ay — 0 
der vor dem Integral stehende Anteil sich zu 
A® errechnet. Bemerkenswerterweise erhilt 
man als Integralwert des diffusiven Transports 
den Wert Null. Es gilt schliesslich : 


Q = 7 wey. 


Das ist gerade die durch den Zulauf zeitlich 
zugefihrte Indikatormenge. 

Aus der Gleichung (3) gewinnt man durch 
geeignete Konzentrationsmessung die Aus- 
breitungskennzahlen Bo, und Bo,.  Allerdings 
dirfte Bo, einfacher durch eine Verdrangungs- 
messung zu bestimmen sein, so dass in diesem 
Fall nur noch Bo, aus der Gleichung (3) zu 
entnehmen ist. 

Unter gewissen Bedingungen wird Gleichung 
(3) von Bo, unabhangig. Entwickelt man namlich 
den Exponenten und den Nenner der Reihenterme 
der Gleichung (3) in Potenzreihen, 


A. 12 ¥ 
4 A,? 


so wird klar, dass die Abhangigkeit von Bo, fir 
einen Einzelterm der Gleichung (3) entfallt, 
wenn in diesen Potenzreihen die zweiten Glieder 
klein gegen Eins sind. Dies (und tiberhaupt die 
Konvergenz der Potenzreihen) kann jedoch nur 
bis zu einem bestimmten A, gelten. Soll nun die 
Unabhangigkeit von Bo, fiir Gleichung (3) als 
Ganzes erfillt sein, so miissen die Terme mit 
A, > A, insgesamt schon klein sein gegen die 
vorigen, was im allgemeinen mit wachsendem 
CA, der Fall ist. 


A,? 
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Kann man 


die genannten 
numerisch realisieren, so vereinfacht sich Gleich- 
ung (3) wie folgt : 


Bedingungen 


R A) 
A A A, 


y=2 


(A 
o ( EXP | (+) 


Bei unserem, in der Einleitung beschriebenen 
Messverfahren wird als mittlere Konzentration 
im Austluss gemessen : 


RA 
y= | 
a 


Mit wachsendem Saulenradius nahert sich der 
Verlauf der Funktion dem  Grenzwert fir 
R A +0, Der Grenziibergang ist nach folgender 
Umformung méglich : 


R * 

R A A * AA * 


exp {| — A,** (Bo,, R® d.L)} AA,* 
[A,* A)) [A,* (RR 
Beim Grenziibergang R A +0 sind nur solche 
von Null verschiedene = A, R méglich, 
fir die A, gegen Unendlich strebt. Bei Grenz- 


tibergingen A + © verhalten sich die Bessel- 
funktionen wie 


J, (A) = sin (A — } 2). 


Mit J,(A,) = 0 folgt J,(A,) = y (2/A, 7) und 
AA, = 2. 


Zur Messung der transversalen effektiven Diffusion in durchstrémten Fiillkérpersiiulen 


Hiermit erhalt man (nach der Definition des 
bestimmten Integrals) 


Derselbe Grenzwert ergibt sich fiir d/A +0. 
Die Berechnung des uneigentlichen Integrals 
erfolgt tiber eine Potenzreihe : 


a 


= (- 1) (n 1)!’ 


(2n + 2)! 
= (-1) 2 nt (mn + 2)'(n + 1) 22 


Folglich : 


(2n + 2)! 
2 Di (- 1) (n + 2)!(n +1) 2 * 


2 
ab exp | Bo, R 


2 
oder wegen | exp - | dé = 
a 


_ (2 +2)! k? 
y 1) 2*"(n + 2)(n Be) 
n-0 


y= 1)" a, (6) 


£ Bo, 


2. Teildarstellung der Tabelle 1. 


Tabelle 1. Querausbreitung bei 
vernachlissigbarer Lingsausbreitung 


| 


0,9000 
1,0000 
1,1000 
1,2000 
1,3000 
1,4000 
1,5000 
1,6000 
1.7000 
1, 8000 
1.9000 
2.0000 


2.5000 
38,0000 
4.0000 
5.0000 


0.3741 
0.4135 
0,476 
0,521 


(2n + 1) 
 2(n+1)(n+2) 


ag = la 


Bo,,- 


Die Funktionsverliufe y wurden aus 
Gleichung (6) fiir den Parameter = 0 
sowie aus Gleichung (5) fiir R A = 11,90 und 
R A = 20/90 berechnet und in Tabelle 1 und 
Diagramm Abb. 2 dargestellt. Mit Hilfe des 
Diagramms findet man die zu einer gemessenen 
mittleren Konzentration gehérende Ausbreitungs- 
kennzahl Bo,. 
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é 
F 
R R R 
=- —_._Bo,,| — = 00000 | — 
0.0000 0.0000 ,0149 0,0494 
0.1000 00244 0.0496 
0.2000 0,0476 0,0476 | 0,0600 
= 0.3000 0.0697 0.0720 
0.4000 0,0908 0.0914 
0.5000 01108 
0.6000 0.1300 0,1300 
OL. 
13 0,7000 0.1483 
0.8000 0.1661 | 
0.1825 
0.1986 
0.2151 
0.2285 
0.2426 | 
0.2560 
0.2690 = | = 
| 0.2813 2 
dL 0,293: Ps 
29333 | 
0 0.3047 
0.3157 
| 
x 
i: 
T 
0 
~ 
| #02222 
W222 
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FORMELZEICHEN 


Zeitlicher Durchfluss des Indikators durch den 

Saulenquerschnitt kg /sec 
Radiale Zylinderkoordinate m 
Radius der Zulauf- und Entnahmerdhre m 


r 


A 


c — Ortlicher Wert der Indikatorkonzentration kg 
Mittlere StrOmungsgeschwindigkeit in 


é — Uber den Querschnitt der Saule gemittelte 
Konzentration kg m3 Achsenrichtung 
rR 


Indikatorkonzentration im Zulauf kg m? x Bo 


L.a 


d — Fir die Packung charakteristische Linge, z.B. 
Teilchengrésse m 
k 


L Abstand zwischen Zutauf- und Entnahmerdhre =m 
A, — Entwicklungsparameter Axiale Zylinderkoordinate 


M, — Longitudinaler Misehungskoeflizient m? see 
M,, — Transversaler Mischungskoeflizient m?® sec d 
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The optimum operating conditions in sequences of stirred tank reactors 


R. Aris 
Department of Chemical Engineering, University of Minnesota, Minneapolis 14, Minnesota 


(Received 8 February 1960) 


Abstract—The problem of choosing the operating conditions, such as temperature and holding 
time, in a sequence of stirred tank reactors so as to maximize the profit from the conversion is 
here solved by a functional equation. This equation, derived by the method of dynamic pro- 
gramming, is very suitable for the repetitive calculations that a modern computing machine 
ean perform. A general algorithm is given both in the case of a cascade of reactors and also 
when there is a by-passing of the original feed to later reactors. 


Résumé —Présentation dune équation fonctionnelle permettant de résoudre le probléme du 
choix des conditions opératoires, telles que température et temps de rétension, dans une succes- 
sion de réacteurs (avec mélange par agitation) pour obtenir un optimum de conversion. Cette 
équation, dérivée par la méthode de programmation dynamique est particuliérement convenable 
pour la répétition des calculs a effectuer avec un computeur moderne. Présentation d'un algorithme 
général convenant au cas dune succession de réacteurs ainsi que lorsque lalimentation primaire 
est partiellement court-circuitée aux réacteurs ultérieurs. 


Zusammenfassung — Das Problem, die Reaktionsbedingungen, wie zum Beispiel die Temperatur 
und die Verweilzeit, in ciner Riihrkesselkaskade so zu wiihlen, dass der Reaktionsfortschritt auf 
einen Maximalwert ansteigt, wird hier durch eine funktionelle Gleichung gelést. Diese Gleichung, 
die mit Hilfe der Methode der dynamischen Programmierung erhalten wurde, ist besonders fiir 
den wiederholenden Rechnungsgang geeignet, den die modernen Rechenmaschinen ausfiihren 
kénnen. Ein allgemeines Zahlensystem wird angegeben, das sowohl fiir den Fall einer Riihrkessel- 
kaskade als auch fiir cine Riihrkesselfolge gilt, bei der cin Teilstrom des urspriinglichen Frischzu- 
gangs den cinzelnen Reaktoren zugefihrt wird. 


Ix aN earlier paper [1] it was shown how the 
problem of choosing the optimum conditions in 
a sequence of stirred tank reactors could be 
solved by dynamic programming. The example 
discussed there was the very elegant one that 
Densicu [2] devised to show the striking improve- 
ment of yield when two stages are held at different 
temperatures. It is the purpose of this paper 
to show how the same methods can be used to 
solve quite general optimum problems for the 
sequence of N stirred tank reactors. The presen- 
tation will be complete in itself, the earlier paper 
serving to give the picturesque detail of a par- 
ticular example. 

More recently Densicu [3] has examined the 
optimum conditions in a sequence of reactors, in 
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which the extent of reaction at any stage is 
controlled by the holding time or the amount of 
reactant added. He achieves the optimum by a 
method of * maximization of rectangles ” drawn 
under a yield curve. It will readily be seen that 
this idea has considerable aflinity with the notion 
of dynamic programming, in fact a simple and 
direct method of calculating his rectangles is to 
use a dynamic programme. 


1. Tue Reacrion System 


We consider a system of m_ simultaneous 
reactions between the chemical species present. 
Since the reactions are independent it is always 
possible to choose a set of m chemical species 
which we will call the key components, and 


R. Arts 


associate each with one of the reactions. From the 
stoicheiometry of the reactions we may express 
the concentrations of all others in terms of these 
and of the reference composition from which 
the reaction mixture is derived. The rates of 
the reactions can be combined to give a function 
R,, namely the rate of change of concentration 
of the i key component due to all the reactions, 
The R, will be functions of the concentrations 
of all reactants and of the temperature, T. By 
the stoicheiometric relations they may be reduced 
to functions of ¢,.¢,...¢, and T and we write 
them T) Ryle, T). The concen- 
trations may be in mole or mass fractions, partial 
pressures or any other units suitable to the par- 
ticular problem. 


%. 


Fic. 1. The reactor system. 


These reactions are to take place in a sequence 
of N reactors such as is shown in Fig. 1. It is 
convenient to number these from the last to the 
fist so that the feed stream enters reactor N 
and the products leave reactor 1. We can quite 
consistently denote the concentrations within and 
leaving reactor » by ¢,, and the feed concentra- 
tions by ¢, y.,. The set of all concentrations 
Cpe t= 1, 2,...m will be denoted by ¢, ; this 
describes the state of the reactant stream leaving 
reactor n and entering reactor (” — 1). In any 
particular stage a choice of the operating variables, 
for example the holding time and temperature, 
must be made and these will affect the way in 
which the reactor transforms the state of the 
stream from ¢,,, to c,. If the flow rate of reac- 
tants is g a balance over the reactor gives 


9 (Cin — = T,,) (1) 


i 


These m equations can be solved simultaneously 
to give c, in terms of ¢,,,. Let this be denoted 
by ¢, F, where 7, is the transforma- 
tion effected by reactor n and depends on the 
choice of the operating variables Vi, and T,,. 


The composition of the product stream can be 
calculated from that of the feed repeated applica- 
tion of this, 


The overall transformation depends on the 2N 
variables 1,2,..N. the case 
of a gas reaction the pressure might be taken as 
an operating variable, but the two chosen here 
will suflice. 

A heat balance over reactor » gives the rate 
at which heat must be added or removed to 
maintain the operating temperature. Thus, if 
c, and p are the specific heat and density of the 
reactants and HH, suitably defined heats of 
reaction, 


2H, R,(c,. T,) — 


qe, p(T, 


where Q,, is the rate of heat removal. Then 


i 


An optimum N-stage policy with respect to 
the state of the feed ey,, is a choice of the 2N 
operating variables that maximizes a certain 
profit function. This profit function will represent 
the gain in value of the products over the raw 
material less the costs of operation. Let v 
be the value of a unit amount of the reaction 
mixture whose composition is Then the 
increase in value of the process stream per unit 
time is 


q {v — 


The operating cost, including the amortization 
of the initial cost, of any reactor is a function of 
its operating conditions, say C(V. In the 
construction of this cost function we must take 
into account the best or current engineering 
practice, which will fix the cost of such items as 
the stirring of a vessel of volume V,, or the heat 
insulation of a reactor operating at 7',. There 
will also be a cost of pumping the reactant steam 
through the reactors and, assuming that it is 
not greatly affected by composition changes, 
this may be taken as a function of g, say C’ (q). 
The profit from an N stage reactor is thus 
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Py lene T3q) = ley) — 
~ 
! 


and this may be calculated by equation (2) 
when V,.... Vy. Ty... Ty and q are chosen. 
Thus when the operating conditions are fixed 
Py becomes a function only of ey,, and Ty,, 
and we may write 


Ty 


The temperature of the feed Ty,, enters this 
expression since it is required to calculate the 
heat removed in reactor N and hence the cost of 
operating this stage. With certain simplified 
profit functions the feed temperature may not be 
important. 

We observe that ¢ enters this equation on rather 
a different basis from the other operating 
variables. This is because it is the same for all 
stages and it is convenient to think of it fixed. 
If we denote by Lp, the profit per unit rate of 
flow without taking the pumping cost into con- 


sideration 


Max P, (5) 


Py =q — 
1 
where 
P, = v(e,) — v — T,) (6) 


then we can determine first the maximum of 
2'p, for any q and afterwards determine the q¢ 
which will maximize Py. 

There will usually be certain physical restric- 
tions on the choice of operating variables, such as 
the requirement that temperatures should not be 
too great. These can be expressed as inequalities 
and we will consider only T,, and V,, subject to 


T, <T,<T* and <V,<V* (7) 


2. Tue Principe or Oprimatiry 


The basic notion of dynamic programming 
has been enshrined by Bettman in the following 
words [4]: “ An optimal policy has the property 
that whatever the initial state and initial decision 
are, the remaining decisions must constitute an 
optimal policy with regard to the state resulting 
from the first decision.” In terms of our reactor 


system we may translate this to read: The 
optimal choice of operating conditions for N 
reactors with respect to the state ey... T'y4, of 
the feed must be such that the conditions in 
reactors (N — 1), (N — 2),...2, constitute 
an optimum (N — 1) stage choice with respect 
to the state cy, Ty. More coloquially, one must 
do the best one can with the output of the first 
reactor if there is to be any hope of achieving the 
best possible result from the complete set of N 
reactors. Moreover, if we know how to choose 
conditions optimally for an (N’ — 1) stage reactor 
we can find the optimal N stage policy by vary- 
ing the conditions in the first reactor and always 
using the optimal (N —1) stage policy with 
respect to the feed from the first reactor to the 
subsequent (N — 1) stages. The maximum over 
all such variations will give the best result for the 
N stage process and in this way we can build up 
the N stage policy a step at a time. At each step 
it is necessary to vary the operating conditions 
in only one stage. Since the amount of work 
required to find a maximum increases exponen- 
tially with the number of simultaneously varied 
quantities the principle of optimality has saved 
a tremendous amount of work and introduced 
a vital organization into the calculation. 
For any particular q let 


N 
= Max (2'p,) (8) 
1 


where the maximization is by correct choice of 
transliteration of the principle of optimality 
gives 


(ewer = Max {py + (ey, T'y)} (9) 


where this time maximization requires only the 
correct choice of Vy and Ty. At any stage we 
can find fy, for 


Ty) = 
Max — C’(g)} (10) 


the maximization being for a correct choice of q. 

Although in analysing an N stage process with 
respect to a particular feed state we are forced 
to consider policies with respect to other feeds 


Chem. Engng. Sci. Vol. 13, No. 2. October, 1960. 


77 


R. Aris 


this extra information is of interest and impor- 
tance. It would be a short-sighted analysis that 
did not consider the effect of variations in the 
feed state on the optimum policy, and though there 
is an optimal value of q it is not always possible 
to use it. 


3. Ture CaLcuLATION AND PRESENTATION OF 


Resvuvrs 
The principle of optimality has given us a 
functional equation to solve. Set out in full the 
system of equations ts 


By = Max {py + (ey. Ty)} 


where 
Py = vley) — — C(Vy. Ty) (12) 
and 


ey Ty (13) 


Fy, is the solution of the equations (1) and in 
calculating the cost function C (Vy, Ty) equation 
(3) is needed to give the cost of heating or cooling 
the reactor. The sequence of equations is started 
by observing that 

and so 


(14) 


The modern digital computer is well suited to re- 
petitive calculations of this type, for it can be pro- 
grammed to do the calculation for given N, print 
out the results, store gy in the place where 
£y-, was before and then go through to same 
eyele of calculation to get gy.,. The machine 
looks for the maximum by direct calculation, 
and inequality restrictions, such as (7), actually 
help by confining the area of search. The chief 
difficulty lies in storing a function of (m + 1) 


variables and this may be prohibitive if m is not 
small. Two devices, both of which are used 
in the earlier paper [1] can be used to cut down 
the work and storage required. Firstly if all 
the reactions are first order the equations (1) 
will be homogeneous in the ¢;. It is then possible 
to work with the ratio of the concentrations to 
any one of them, say ¢,,, and this reduces the 
number of variables by one. Secondly, it may 
be that the temperature enters the rate ex- 
pressions R, only in the form of a rate constant 
k,exp(— E, RT). In this case if we are not 
concerned with costs dependent on temperature 
it is convenient to take one of these rate constants 
as an independent variable. All the other con- 
stants can be expressed as powers of this one and 
the maximizations with respect to 7’ may be 
performed by varying this rate constant. This is 
perfectly satisfactory since exp(— E RT) is a 
monotonic function of 7. 

The results of each step of the calculation may 
be tabulated as shown below : 
These tables are prepared one after the other by 
solving the equations (11), (12) and (13), Given 
any feed cy... Ty., we can find the complete 
optimum policy as follows. First by entering the 
table for n= N at ey, Ty, we have the 
optimum profit function gy (¢y.,. Ty.) and the 
operating conditions, Vy, Ty, for the first stage. 
The final columns give the state of the stream 
produced under these optimum conditions, Since 
the principle of optimality states that the remain- 
ing (N — 1) stages must use optimal policy we 
can find the operating conditions in the second 
reactor by entering the table for n — N — 1 at 
cy. Ty. This will give Vy, Ty, and the 
resulting state cy ,, Ty , which is the feed to the 
remaining (N — 2) stages. Passing through the 


Optimum profit 
function 
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tables in this way the complete optimum design 
may be rapidly extracted. 


4. Tue Norion or toe LaGrance 

In the general profit function described above 
a definite cost was associated with the volume of 
the reactor. It might happen however that the 
interesting problem is to determine the maximum 
yield of a particular component with little or no 
regard to the costs of operation. In some cases 
a physical restriction will still have to be imposed. 
For example with a reversible reaction the greatest 
yield will be the equilibrium yield at the most 
favourable temperature but it will require an 
infinite volume of reactor to achieve this. We 
are thus interested in the optimum yield for a 
fixed total reactor volume V. This can be 
expressed by explicitly letting the optimum yield 
be a function of V; thus if ¢ is the concentration 
of the product we let 


Sy V) = Max fey 


the maximization being over choice’ of 
T,, T,,... Ty and V,, Vy... Vy subject to 


The principle of optimality then gives 


(exe V) = 
Max {ey — + fy-alew, V — (17) 


the maximization being over choice of Ty and Vy 
subject to Vy < V. Here, however, the function 
Sy is a function of two arguments and the storage 
and interpolation of such a function is much 
more troublesome than that of a single variable. 

Another approach to this problem is to con- 
struct an artificial cost function associated with 
the volume and to maximize this slightly more 
general profit function. Thus in this case we 
would consider the profit function 


(15) 


(16) 


N N 
where 


and yz is any positive constant. Now if 


Pa =n — — 
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we have no need to make gy a function of J V, 
for the maximum will occur at a finite total 


N 
volume. In fact now V = Z V,, is a function of 
1 


#, for the equation for gy, namely 


= Max {py + (€n)} (20) 


can be solved for any N, » and ey,, and the 
total volume calculated. For a given feed state 
V is such a function of » and N as is shown in 
Fig. 2. It is now possible to choose the value of 
p# so that the total volume will be as required. 
pu is called the Lagrange multiplier and is in fact 
a cost associated with the volume so that the 
larger « is the smaller will be the optimum volume. 
This way of putting the problem is often of 
great advantage in saving computer storage and 
time. 


Fic. 2. Dependence of the total volume on the 
Lagrange multiplier ,. 


The Lagrange multiplier can also be adapted to 
other situations. If, for example, we seek the 
policy for an N stage process which gives the 


N 
greatest space-time yield (ec, — ey,,)/2 V,, it is 


computationally easier to solve the problem 
defined by (18), (19) and (20) and afterwards 
find the value of « which maximizes the space— 
time yield. The same technique can be employed 
if we seek to maximize the fraction of one of the 
reactants which is converted into a particular 
product. The value function v(c) can also be 
used to include specification requirements on the 
product. Thus if the concentration of a particular 
component is required to lie in a certain interval 
v (c) can be put equal to zero outside this interval. 


Chem. Engng. Sci. Vol. 18, No. 2. October, 1960, 


R. Ants 


5. Pouicy wirn By-passinc or Freep 


As an extension of the general programme it is 
interesting to see how to find the optimum policy 
with by-passing of feed to later stages of the 
process. The system is as shown in Fig. 3, at stage 
n we have not only the product of the previous 


Cores 


Fic. 3. The reactor system with by-passing. 


stage (n + 1) coming in with flow rate q,.,, but 
also added, fresh feed of fixed composition cy at 
a rate q,’. The mass and heat balances now read 


= Un+1 (21) 


Qn Cin — Cintas — Cie = VR (e,. T,,) (22 
V, Rew — (23) 
1 


where h; = H, ¢,p and q,@,* = @, The 
operating conditions at each stage are now not 
only V,, and but also q,’. 

We may again write a general profit function 
in the same form as (4) above : 


Py Trai Co To; qi 
=, 0 — ? — 


N N 
v — Vie Ty) — (24) 


where C” (q,) is the estimate of the cost of pumping 
in terms of the final flow rate. As before q, 
enters in rather a different fashion from the other 
variables since all else depends only on the ratio 
of the flow rates. We shall again adopt the 
technique of seeking first the maximum of 


N 
{Px + aq Lp, and afterwards choos- 


ing q, to maximize Py. At the expense of a 
little more complication however it would not be 
difficult to include the individual pumping costs 
of all side streams. 


Let = An =n = An — and 
Pu = AgVlCn) — Ansa (Cass) — 
— A,’ v (eg) — {C(V,,. T,) 
= A, {v (e,) — v — Ansa 
Then if 
Ty. 13 Co To) = Max Sp, (26) 


(25) 


the maximization is over suitable choice of 
V,,... Vy, Ty... Ty and Ag, .. . subject to 
S Ay S AQ SAY = 1. We should not 
at this stage identify the feed composition 
€xy., With eg for we have to set up a functional 
equation for gy as a function of its first two 
arguments. At the final stage we can identify 
Cys, and cy and then of course Ay,,, loses its 
meaning. The presentation also allows for the 
by-pass feeds being of different composition from 
the original feed: it would be possible to have 
arbitrary compositions for all these feeds but 
it would involve carrying along many more 
variables. When the meaning is quite clear we 
can write gy as a function of its first two arguments 
only. 


If we write — Vg, as the holding time 
of stage n, then equations (22) and (23) are 


An (Cin Cio) Anta — Cio) t, Ri(c,, T,,) (27) 


t, Lh, R; (Cys T,,) (28) 

Applying the principle of optimality to gy we 

have 


Tye) Max {Pw + By-a (en: Ty)} (29) 


where the maximization is over suitable choice of 
Ty. ty and Ay,,. Ty and ty may be restricted by 
such inequalities as (7) and Ay, is less than or 
equal to Ay which is known from the calculation 
of gy, The sequence of equations is started 
with 

£1 = Max {p,}, A, = 1. 


The case where the only heat supplied is in 
initially heating the feed to reactor N, and a 


Cw Coe C2 T Gq 
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falling temperature sequence is obtained by by- 
passing cold feed, is of particular simplicity and 
interest ; it will be considered in detail in a later 
paper [5]. Sufficient has been said here to show 
the suitability of dynamic programming for the 
whole range of optimum problems with sequences 
of stirred tanks. 


Acknowledgement—1 indebted Professor 
K. G. Deneicu for his interest in this topic and for 
showing me his contribution to the Amsterdam con- 
ference whilst it was still in manuscript. 


NOTATION 


C(V,. T,) cost of operating reactor as a fune- 
tion of its volume and temperature 
C’ (q) = cost of pumping at rate q 
¢; = concentration of it” key species 
Cp» ~ Specitic heat of reaction mixture 
In len. maximum profit from N reactors with 


feed state Cn. r Ty. 1 


in ley. 
N 


Py 


‘ = by-pass flow rate to n' 


Max p,, 
- heat of it” reaction 
= number of reactors 


= 4 (ey) Bley. 
N 


v(¢,) — 0 — {C(V,, T,) 
- rate of heat removal from n'" reactor 
flow rate through all reactors 
(sections 1—4) 
flow from reactor n 
reactor 
= rate of change of c; due to all reactions 
- temperature of reactor n 
V4. holding time of reactor n 
= volume of reactor n 
value of unit volume of reaction 
mixture of composition c 
= Vy 
= An Anes 


p = density of reaction mixture 
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Abstract—A process is described for the separation of binary liquid mixtures by a method which 
is fundamentally the same as extractive distillation but uses a finely divided solid as a carrier for 
the non-volatile solvent. A non-volatile solvent absorbed on the surface of crushed kieselguhr 
flows down a column countercurrent to a stream of air, The mixture to be separated is injected 
continuously as a liquid feed into the mid-point of the column. The flows of gas and solid are 
adjusted such that one component is carried upwards in the gas stream and the other component 
downwards absorbed in the non-volatile liquid. Benzene and cyclohexane have been successfully 
separated at a rate of 30cm hr in a 1 in, diameter column 4 ft long when using polyoxyethylene 
400 di-ricinoleate as the non-volatile solvent. The H.T.U. (height of transfer unit) for this system 
and column packing was found to be approximately 2 em, 


Résumé Les auteurs décrivent un procédé de séparation de meélanges liquides binaires par 
une méthode semblable a celle de la distillation extractive mais utilisant un solide finement divise 
comme support du solvant non volatil, Un solvent non volatil absorbé par du kieselguhr concassé 
s'écoule du haut vers le bas d'une colonne, a contre-courant d'un écoulement dair. Le mélange 
a séparer est injecté de fagon continue dans le milieu de la colonne. Les débits du gaz ct du solide 
sont établis de telle sorte qu'un composant est entrainé vers le haut par le courant gazeux tandis 
que, Tautre absorbé par le liquide non volatil, s‘écoule vers le bas. Une bonne séparation de 
benzéne et de cyclohexane a été faite, a une vitesse de 30.em* h dans une colonne de 1 pouce 
de diamétre et 4 pieds de haut en utilisant comme solvant non volatil du polyoxyethyléne 400 di- 
ricinoleate. La hauteur dunité de transfert (HITU) pour ce systéme et le garnissage de cette 
colonne est de 2m environ. 


Zusammenfassung —Es wird cin Prozess zur Trennung binirer Flissigkeitsgemische bese- 
hrieben. Diese Trennung wird durch eine Methode erreicht, die grundsiitzlich dasselbe wie eine 
extraktive Destillation darstellt, bei der aber cin fein verteilter Feststoff als Traiger fiir ein 
nichtfliichtiges Losungsmittel dient. Dieses nichtflichtige Losungsmittel, das auf der Oberfliche 
von gemahlener Kieselgur adsorbiert ist, fliesst in ciner Séule entgegen cinem Luftstrom nach 
unten. Das zu trennende Gemisch wird kontinuierlich als Flissigkeit im Mittelpunkt der Siiule 
eingefiihrt. Der Gas- und Flissigkeitsstrom wird so eingestellt, dass cine Komponente im 
Gasstrom nach oben und die andere in der nichtflichtigen Fliissigkeit absorbiert nach unten 
gefihrt wird. Benzol und Cyclohexan wurden erfolgreich getrennt bei einer Zuflussgeschwindigkeit 
von 30 em? hin einer 120 em langen Saule von 2,5 em Durchmesser. Dabei wurde Polyoxvaethylen 
(400)-di-ricinoleat als nichtflichtiges Losungsmittel verwendet. Es wurde gefunden, dass die 
HTU (Hohe einer Chertragungseinheit) fir dieses System bei der angegebenen Siulenfiillung 
etwa betriigt. 


The different constituents are eluted from the 
column by the moving gas stream. A detector 
placed at the outlet from the column indicates 
the time of elution and also the concentration of 
the component thus enabling qualitative and 
quantitative analyses to be made. 


1. INTRODUCTION 


Mvcu interest has been aroused in recent years 
in the use of gas—liquid chromatography (G—L_C) 
as a means of analysing mixtures of volatile liquids 
and gases. In this technique the mixture to be 
analysed is vaporized into a stream of inert gas 


which passes through a column containing a non- 
volatile liquid absorbed on the surface of a finely 
divided solid. Due to the differences in volatility 
of the solute and or physical interaction in the 
non-volatile liquid phase, separation into pure 
constituents occurs as bands within the columin. 


The use of large seale G-L-C columns up to 
% in. diameter to increase the quantity to be 
separated at any one time has been fully in- 
vestigated [1, 2, 3, 4]. A semi-continuous unit has 
been produced commercially in which the next 
batch of material is injected into the column as 
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the last component from the previous batch is 
eluted [5]. 

The purpose of this work has been to study a 
countercurrent analogue of G-L-C. That is to 
allow a solid containing an absorbed solvent to 
move countercurrently to an inert gas stream and 
to inject a liquid feed continuously into the mid- 
point of the column. By using this technique it 
was hoped to achieve diflicult separations of 
volatile liquids and gases on a much greater scale 
than has been possible with G-L-C columns. 

This countercurrent analogue of G_L-C can be 
considered to be extractive distillation [6] since 
the solid acts merely as a carrier and takes no 
part in the diffusional processes occuring. A 
continuous process operating in the above manner 
and using a vapour feed has been used to separate 
aromatic hydrocarbons from coal gas [7]. The 
countercurrent analogue of gas adsorption chrom- 
atography, hypersorption, [8, 9] has been practised 
on a commercial scale for just over a decade. 


2. Turory 


This method, like other processes for separating 
different compounds is dependent upon the 
components distributing themselves differently 
between two phases. In this case the phases are 
air and a non-volatile solvent absorbed on the 
surface of suitable solid such as kieselguhr. The 
solid acts as a carrier and takes no part in the 
separation. 

Each phase has a preferential selectivity for 
one or other of the solute components. Movement 
of two phases countercurrently will produce in 


each phase a relatively high concentration of one 
solute and a relatively low concentration of the 
other, so that given a sullicient height of column 
a substantially complete separation is possible. 

The suitability of a solvent for a given separa- 
tion together with conditions of operation may 
be seen from a chromatogram of the binary 
mixture (AB) obtained on a gas-liquid chromato- 
graphy apparatus (Fig. 1). 

The linear velocity of component A through the 
column is z R, where zis the length of the column 
and R, the retention time of component A. 
Similarly the velocity of B is z Ry. Thus if the 
column packing is moved countercurrently to the 
carrier gas stream at a velocity greater than = R, 
but less than z Ry, component B will be carried 
in the direction of the packing absorbed in the 
liquid phese. Component A will be eluted in the 
carrier gas. 

It is more convenient to express the conditions 
of separation in terms of the volume flow of solvent 
and carrier gas rather than linear velocities. 

The partition coeflicient, k, for a particular 
component may be detined as 


Weight of component per unit volume of solvent phase 
Weight of component per unit volume of carrier gas phase 
Let the partition coetlicients of components A 
and B at the temperature of the rectification 
section of the column be k} and kf respectively 


G 


= 
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(see Fig. 2). If component B has the longer reten- 
tion time in the stationary column, then kf will 
be greater than k*. 

Let the volume flow rates of solvent and 
earrier gas in the rectification section be L and 
G* respectively. These may be considered 
constant when the temperature of the column is 
kept constant, and the pressure drop across the 
column is small. 

For component A to move in the direction of 
the gas phase G" L must be greater than k ¥. 
Similarly for component B to move in the direc- 
tion of the liquid phase, G" L must be less than k,*. 

The condition for separation is 


R 
Thus there will be a range of possible operating 
values of Gp L over which separation may be 


kif < 


obtained. In the ideal case complete separation 
would be achieved over the whole range of values 
of Gy L from kf to kj The actual operating 
values may only approach the respective & values, 
the degree of approach being dependent upon the 
column size and its efliciency. 

For component B to be stripped from the liquid 
phase, the ratio of carrier gas flow to solvent tlow 
in the stripping section for the ideal case is given 
by 


where k,> is the partition coetlicie.:. of component 
B at the temperature of operation of the stripping 
section, and G® is the carrier gas flow rate in that 
section. The operating value of G® L must be 
much greater than kj? in order to recover all of 
component B from the liquid phase. The actual 
value will again depend on the column size and 
its elliciency. 


3. EXPERIMENTAL 
Materials used 


One of the advantages of this process over 
simple distillation is that separations can be 
achieved by differences in polarity as well as 
boiling point. Polar molecules are those in which 
the charge density is not uniform over the surface 


of the molecule but is concentrated at one or two 
points. The method can therefore appropriately 
be demonstrated by separating two liquids with 
approximately the same boiling point, but of 
different polarity. With this in mind preliminary 
work has been carried out using a polyglycol 
derivative (polyoxyethylene 400° di-ricinoleate) 
absorbed on 10 20 mesh kieselguhr to separate 
benzene (b.p. 80-1-C) and cyclohexane (b.p. 
80-7 


Equilibrium data determinations 


To permit the estimation of the relative flow 
rates of packing and carrier gas the partition 
coellicients of benzene and cyclohexane between 
the solvent and air at various temperatures were 
determined, Measurements of retention times 
were made on a standard G-L-C unit, using a 
4mm diameter column, 5 ft in length, packed 
with 40-60 mesh C22 fire-brick containing 25 
per cent by weight of polyoxyethylene 400 
di-ricinoleate, Air was used as the carrier gas, 
with a katharometer as the detecting element. 
The partition coeflicients were calculated from the 
following equation derived by James and Martin 
{10}. 

A, | 
Ay 


(P, 

— | 

cross-sectional area occupied by the 

mobile phase 

Ay — cross-sectional area occupied by the 
stationary phase 

Vi, — retention volume to the centre of the 
Zone 

= — length of column 
P, — pressure at column inlet 
P, = pressure at column outlet 


where A, : 


Values of the partition coeflicients at various 
temperatures are shown in Table 1. It was found 
that the partition coetlicients could be expressed 
by the following equations : 

1440 


For benzene log,)k = 1-975 


For cyclohexane log yok = — 1719 


where k = the partition coefficient, @ = “Kelvin. 
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Table 
Ratio 
Temperature values for K values for 

benzene cyclohexane benzene 
cyclohexane 

oes 208 209 

450 159 2-70 

128 a7 224 

75 : 2.403 


General description of apparatus 


The flow diagram of the apparatus is shown 
in Fig. 3. A sealed reservior (A) containing the 
packing was placed at the top of the 1 in. diameter 
column, the flow being controlled by means of a 
feed table (B) at the bottom of the column. Air 
was supplied to the column below the stripping 
section (C) and was drawn through under vacuum. 
The air was dried, preheated in (D) to aid the strip- 
ping and supplied at a pressure slightly above 
atmospheric, This ensured a flow of air out 
from the bottom of the column so preventing any 
moisture entering the system. 


Inyector 
pont 


L 
Fie, 3. Flow diagram of apparatus 
T Temperature measurement 
P Pressure temperature 
F Flow measurement 


The mixture to be separated was injected at 
the mid-point of the rectification section (E). 
The more strongly absorbed component was 
carried down in the liquid phase and was removed 
from the column by an excess of air in the 
stripping section (C). The less strongly absorbed 
component was carried upwards in the gas phase. 
Both streams passed through cold traps to con- 
dense out the separated products. 


Column construction 

The glass column was divided into four sections, 
the rectification section (FE) being 4 ft in length 
and the stripping section (C) 2 ft. The section 
below the stripper (G) which allowed the packing 
to cool before leaving the column was 18 in. in 
length. The stripping section was enclosed in an 
air jacket contained by the column and a con- 
centric glass tube, 2in. in diameter, on which 
several turns of resistance wire were wound. A 
further concentric tube, 3 in. in diameter, served 
to eliminate the effect of draughts. 


Flow of packing 

The support used was Johns-Manville C22 
firebrick ground approximately to 10-20 B.S. 
mesh and carrying 249 per cent wow of liquid 
phase. The sieve analysis at the commencement 
of the runs is shown in Table 2. 

The packing flowed under gravity from a sealed 
reservoir which contained supply lasting 
4-6 hr, depending on the flow rate. The rate of 
discharge was independent of the head unless the 
reservoir was nearly empty. To avoid irregular 
flow of the packing it was necessary to vibrate 
the column. This was achieved by means of an 
electric bell impinging on the column just below 
the stripping section. 

The packing discharged through an orifice at 
the bottom of the column on to a rotating table, 
which was driven through a gear system by a 
variable speed motor. A range of speeds from 
1-4 rev, min was possible. The packing was re- 
moved from the table by a scraper. The position 
of the table relative to the column outlet could 
be varied, so increasing or decreasing the rate of 
removal of the packing. 

Details of the orifice at the bottom of the 
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column are shown in Fig. 4. The flow rates were 
found to be too high when the width of the solid 
band removed was greater than 0-5 in. and the 
conical section was introduced to reduce the 
diameter of the column. The height of the orifice 
above the table could be varied to control the 
depth of packing removed and orifices ranging 
from 0-25-0-5 in. were used to vary the width of 
the solid band. 


Fig. 4. 


Thus the flow rate could be controlled or 


varied in the following ways: 


(a) Speed of the feed table. 

(b) Position of the feed table relative to the 
column outlet. 

(c) Height of the orifice above the feed table. 

(d) Size of the orifice. 


Air supply flow, temperature and pressure 

measurement 

Silica gel was used to remove any moisture in 
the entering air which was preheated to aid the 
stripping process. 

The flow of air to the column was indicated 
by a rotameter between the drier and heater. The 
temperature and pressure were measured at the 


column inlet. The flows of air from the column 


were measured by capillary flowmeters. The 
temperature and pressure of both streams were 
measured at the column outlets. 
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Injection and collection 


The mixture to be separated was injected by 
means of a micro-pump through a hole, , in. 
in diameter, in a brass flange which was inserted 
between the two sections of the rectification 
column, A small heater around the inlet feed pipe 
ensured vaporization of the liquid into the air 
stream. The maximum capacity of the micro- 
pump was 40 hr. 

The products were recovered in cold traps 
which could be by-passed during the initial stages 
of a run, each outlet stream from the column 
passing through two such traps. The first trap 
contained solid carbon dioxide, the second solid 
carbon dioxide in acetone. The first trap at a 
temperature of — 20°C condensed out most of 
the liquid but a second trap at 80 -C was 
was necessary to recover 95 per cent of the material 
fed into the column, the remainder passing 
through the traps to the vacuum pump. 


Table 2 


Sieve analysis of brick after 
Sieve analysis of starting being recycled approximately 


material eight times 
B.S. mesh Wt. B.S. mesh 
Ol4 
s 10 32-00 18-51 
12 13-23 956 
12-14 22-47 12-14 19-96 
14-16 2140 14-16 32-28 
16018 160-18 17-27 
Is 25 18-25 


Experimental procedure 


The air jacket surrounding the stripping column 
was allowed to attain the required temperature, 
with air flowing through the column. The feed 
table was then started and the gas flows adjusted 
to the required values, The pressure at the inlet 
to the stripper was maintained at approxinately 
5em water-gauge. Once the flows were steady 
injection was commenced, The apparatus was 
run for about 30 min to reach steady conditions 


after which the outlet streams from the column 
\ 
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were diverted through the cold traps. Since 
sample volumes cf 10 were sullicient for 
analytical purposes the length of each run de- 
pended on the injection rate. The solid support 
was collected throughout this period and the 
average flow rate determined. 


Analysis 

The samples were analysed by gas-liquid 
chromatography, using the column described 
previously, but with nitrogen as the carrier gas. 
By increasing the sample size concentrations 
down to 0-1 per cent could be detected. 


4. Resvuvrs ann Discussion 


Complete separation of benzene and cyclohexane 
has been achieved over a wide range of operating 
conditions, the results being summarized in Fig. 
5. The lines 4B and CD show the limiting values 
of air and solid tlows for complete separation. 
Complete data for three of the runs are shown 
in Table 3. The injection rate was maintained at 
30 em® hr throughout this series of experiments ; 
the air flow and temperatures through the 
stripping section also remained constant. 

At an injection rate of 30cem* hr and an air 
flow rate of 501 hr the air was saturated with 
vapour around the feed point, and the packing 
appeared moist for some distance below the in- 


~ complete separation 
© cyclohexane in bottom product 
[) benzene in top product. 
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jection point, though separation was accomplished 
at lower air rates. The flow rates of air and solid 
obtainable with the present design limit the 
throughput to 40 em® hr. At higher air and 
solid rates, complete separation occurred over a 
range of GL values from 300 to 750. As saturation 
was approached this range was reduced to 350 to 
700, The limiting values of G L for separation 
were predicted as 300 and 900 at 20°C as can be 
seen from Table 1. Benzene appeared in the top 


Table 3 


2 


Run number 1 


|" 


Air flow in rectification section (1/ be) 115-7 


101-3 


Air tlow from stripping section nr) | 109-9 106-3 


a hr) 


Air tlow to column 257 257 


Outlet pressure ; rectification section 
(mm He) 744 739-5 


Outlet pressure THES 746-3 | 


: stripping section 


764-9 | 764-8 | 764-2 


Inlet pressure : stripping section 
(mm Hg) 


Solid flow ig hr) 


Injection rate (om? 


cc 91 
_ 
97 


v2 


Temperatures 7, 


Overhead product : 
(% vol.) 


999 99-9 | 


99-2, 99-9 


Bottom (% vol.) | 


product when GL was approximately 700, 
showing that the number of theoretical stages 
above the injection point was insuflicient to 
remove all the benzene from the gas phase. 
Below the injection point the rectification section 
followed the predicted value of GL, ite. 
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cyclohexane appeared in the bottom product when 
G L was approximately 300, The apparent higher 
efliciency in this section was due to the fact that 
the air entered the rectification section at a 
temperature of about 60°C. The partition co- 
eflicient for cyclohexane at this temperature was 
87, and the condition for cyclohexane to move in 
the direction of the gas phase was that the G L 
be greater than 87. The operating range over 
which separation was achieved was thus increased 
by the temperature gradient existing at the 
bottom of the rectification section. 

It is apparent that a rectification section of 4 ft 
in length is capable of separating a system having 
a partition coeflicient ratio less than that of 
benzene cyclohexane. 


Evaluation of column efficiency 


By definition for dilute systems * the number 
of over-all gas phase transfer units.” Nog. with 
reference to the heavy component benzene, is 
given by equation (1) 

Pa 
Noo = 
For a particular absorption column of packed 
height Z the number of transfer units and the 
height of a transfer unit are related 


Z = Hog * Noe (2) 


If the equilibrium line is linear Nog can be 
calculated from the terminal concentrations 


> (Wy Ye) 


No (3) 
(Y Ye) LM. 


For the experimental conditions of this in- 
vestigation the resistance to mass transfer resides 
in the gas phase due to the high solubility of the 
solutes in the solvent; any resistance in the 
liquid phase may be neglected. 


Under these conditions 


Hog = He 


Nog = Ne 


To illustrate the method a set of conditions has 
been chosen from the line AB on Fig. 5 which 
represents the upper limit of air flow for a given 


solid flow, i.e. benzene was beginning to appear in 
the top product and the column above the feed 
could be assumed to be working completely. A 
temperature of 20°C and an injection rate of 
30em* hr have been assumed similar to the 
experiments recorded in Fig. 5. Shown plotted 
on Fig. 6 are the equilibrium line for benzene at 
20°C and the operating line for benzene in the 
section of the column above the feed point, at 
column conditions indicated in Table 4. The 
purity of the cyclohexane at the top of the column 
has been taken as 99-9 per cent. 


Slope 1/889 


Concentration in gos phase, 


Air flow rate 


Solvent flow rate 


105-2 


(mi/hr) 143-6 


(ghr) 


cyclohexane (g hr) 11-57 


Injection rate : benzene 


(gem) 56 x 10°5 


cyclohexane at feed point (g/em*) 78 x lo 


benzene in top product (g cm?) | 1x 107 


cyclohexane in top product (g em*)| 1-1 10-4 
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Because of the high solvent concentrations and 
isothermal conditions existing above the feed 
point it was reasonable to assume that the 
equilibrium relationship was linear. Also the 
operating line for the benzene is linear when con- 
centration are expressed in the units shown. 

The value of Ng, between the top of the column 
and the feed point by equation (3) was found to 
be 29-9 for the conditions given in Table 4, which 
corresponds to an Hg value of 22cm. This 
method was not applied to the section of the 
column below the feed point because of the 
unknown temperature gradients there. 
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NOMENCLATURE 


Ay, = cross-sectional area occupied by the mobile phase 
en? 


Ag = cross-sectional area occupied by the stationary 
phase (em?) 
R = retention time to centre of zone (min) 
=z = length of column (cm) 
Vp — retention volume to centre of zone (ml) 
Pj ~ pressure at column inlet (mm Hg) 
Po ~ pressure at column outlet (mm Hg) 
k = partition coefficient 
L = solvent flow (1/hr) 
G = air flow (1 hr) 
H¢, = height of a gas film transfer unit (em) 
Hog, = height of an over-all gas phase transfer unit (em) 
Ng = number of gas film transfer units 
Nog = number of over-all gas phase transfer units 
y concentration of component in gas phase (g 
Subscripts 
A = component A 
B = component B 
1 =concentrated end of countercurrent system 
2 = dilute end of countercurrent system 
e = equilibrium value 
L.M. = logarithmic mean 
Superscripts 
R = rectification section 
S = stripping section 
F = feed point 
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Letters to the Editors 


Effects of sonic vibrations on the rate of absorption of gases from bubble beds 


Ir 1s well known that bubbles produced at jets and in 
boiling liquids emit characteristic sounds. Mixnacer [1] 
first analysed the sounds of bubbles produced at jets 
under various liquids and concluded that most of the 
sound was produced by radial pulsation of the bubble. 
Shortly afterwards Swirn [2] considered the reverse 
effect the vibration of bubbles produced by sound 
generated in the liquid — and concluded that the bubbles 
should show a high degree of resonance since the damping 
coeflicient was : ‘nall. He also suggested that the intense 
local strains in jhe vicinity of a pulsating bubble might 
account for the destructive effects of sound on vegetation 
and other living matter. In this connexion AcKERMAN 
{3} has shown that acoustic vibrations in the frequency 
range 0-4-3500 ke/s can destroy Paramecia and Mammalian 
erythrocytes. The cavitation which destroyed the cells 
was ascribed to the presence of extremely small bubbles. 
Cavitation induced by sonics and ultrasonics in liquids 
has also been discussed by Rust [4], Karusrin [5] and 
Meyer [6). In view of the above work it seemed likely 
that bubble vibration phenomena could be used to improve 
mass and heat transfer in operations such as gas absorption 
and boiling. As part of an initial study in this field we 
have studied the effect of sonic vibrations on the absorp- 
tion of carbon dioxide from bubbles in water. 

The apparatus is a modification of that described by 
Hovcuron et al. [7] and is shown in Fig. 1. The gas 
and liquid phases continuously enter and leave the top 
of a 2 ft length of 3 in. diameter glass pipe. The bottom 
of the column contains a closely fitting polystyrene 
piston sealed by a rubber diaphragm. The piston was 
vibrated by an electrodynamically driven vibrator at 
frequencies in the range 20-2000 ¢ ss. The amplitude of the 
vibrations could be measured at each frequency by means 
of a microscope. A plot of amplitude against frequency is 
shown in Fig. 2. In order to avoid the effects of inert 
gases on the rate of absorption pure carbon dioxide was 
used. The gas was distributed in the liquid by a sintered 
glass candle of medium porosity. The rate of absorption 
was calculated in terms of K,;a by measuring the inlet 
and outlet gas and liquid rates and also by analysing the 
liquid phase for dissolved carbon dioxide. 

Fig. 2 shows the effects of frequency and amplitude 
on the rate of absorption of carbon dioxide in water as 
measured by the K, a. The measurements were made at 
constant temperature, pressure, bed height, gas flow rate 
and liquid flow rate. It is evident from the data that the 
rate of absorption can be increased by as much as 50 per 
cent when certain vibration frequencies are used to excite 
the bubble bed. The greatest effect is at about 75 ¢ 's. 


Carbon 
donde 


Elec tro- dynormc 
vibrator 


Fic. 1. Equipment. 


The fact that the absorption rate shows peaks at certain 
frequencies together with the relatively small amounts of 
power introduced into the liquid (less than 4W at all 
frequencies) indicates that the phenomena are due to 
resonance effects. The plot of amplitude vs. frequency for 
the case where no bubbles are present shows that peaks 
in amplitude cannot be the cause of the peaks in absorp- 
tion rate. In most cases the peaks in absorption rate 
occur when the power input to the liquid is decreasing or 
relatively constant. Resonance effects are also indicated 
by the fact that the amplitude of the forcing vibration is 
somewhat dependent upon whether bubbles are present 
in liquid or not, as shown in Fig. 2. This effect is also 
quite audibly noticeable. When bubbles are introduced 
while the column of liquid is being vibrated, the presence 
of bubbles seems to strikingly increase or diminish 
(depending on the frequency) the volume of sound heard. 
Whether the effect as measured by K, a is due to a change 
in film coefficient for the individual bubbles or due to a 
change in area for absorption is not clear at the moment. 
However, visual observations on the effects of sound on 
single bubbles show that the sound can greatly affect 
their appearance. For instance, at low frequencies 
(20-60 ¢ s) the effect of sound appears to elongate the 
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Fic. 2. Effect of frequency and amplitude on K,a for CO, in water. 


Pressure = 730 mm Hy. 
velocity — 0-050 ft see. 


bubbles in the direction of vibration while at higher 
frequencies (800-1200 ¢ s) the effect of the sound is to 
considerably disturb the bubble boundary to such an 
extent that the bubble takes on a frothy appearance. 
Intermediate frequencies (40-800 © s) seem to have little 
effect on the appearance of single bubbles but greatly 
affect the appearance of the bubble bed, especially at 
the intermediate peaks shown in Fig. 2. For instance, 
at the peak of 75 ¢/s the visual effect of the sound is to 
produce turbulence and coalescence in the bubble bed. 
At the surface of the liquid the bubbles burst with some 
violence. Possibly these effects are caused by the presence 
of pressure nodes and antinodes as discussed by GOLDMAN 
and RinGo [8] and also Biake [9]. It is believed that the 
peak at about 840 ¢ s shown in Fig. 2 might correspond 
to the pulsation frequency of Minnarrr [1]. Since the 
pulsation frequency depends upon the diameter of the 
bubbles the average diameter of the bubbles in the bed 
ean be calculated from Mixxarrr’s formula below 


656 

f D 
where f is the frequency in ¢/s and D is the bubble»dia- 
meter in cm. The diameter corresponding to 840 c¢ s 
would be 7-8 mm. This is a much larger diameter than that 
usually found under similar conditions without the presence 
of sound [7] thus indicating that the sound may have 
produced coalescence of bubbles in the bed, thereby 


Temperature = 15°C. 
Water velocity 


Inlet superficial gas 


0-040 ft sec. Bed height — 1-0 ft. 


reducing the area for mass transfer. If this is the case the 
film coefficient must have increased considerably in 
order to produce an over-all increase in K,; a. The effect 
of sound on the coalescence of bubbles is discussed by 
Biake [9], who points out that two bubbles would be 
attracted to each other if their oscillations were in phase 
and repulsed if their oscillations were out of phase. It is 
possible that excessive coalescence might have caused 
K,,4 to decrease below the K, a for no vibration shown at 
240 ¢ sin Fig. 2. The presence of peaks at low frequencies 
indicates that bubbles can vibrate in modes other than 
pulsation. However, the peaks at 145 and 280 c/s may be 
overtones of the peak at 75 .¢/s since they are of smaller 
intensity. Recent work by Kosre.ia [10] and SKLENCAR 
{11] using a wave-analyser has confirmed the existence 
of the pulsation frequency of Minnarrr, but has shown 
that frequencies other than the pulsation frequency can 
be emitted from single bubbles formed at a jet. Although 
these workers did not make measurements below 200 ¢ ‘s 
they found other frequencies that could not be attributed 
to pulsation, particularly in the range 400-600 
Although no attempt was made to attribute these fre- 
quencies to definite modes of vibration, they may account 
for the generally high level for the absorption rate in the 
400-600 © s region shown in Fig. 2. 


K. L. Harsaum 
G. Hovcnron 
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Book Reviews 


F. C. and C, E. Drypen: Chemical 
Engineering Plant Design. (Fourth Edition — revised). 
McGraw-Hill, New York 1959. 534 pp. 80s., $12.00. 


Since the inital publication of Chemical Engineering 
Plant Design by F.C. Vicepranpr in 1934, this book has 
been widely used through three editions as both a text 
and a reference work for plant design in the process 
industries. 

The fourth edition, under the co-authorship of F. C, 
Virpranpr and Drypen, presents a considerable 
departure, both in content and arrangement, from the 
previous editions. Information is presented in approx- 
imately the same chronological order as it is required for 
the actual design of a chemical plant, beginning with a 
consideration of development of the project, and proceeding 
through process design, equipment and materials selection, 
plant layout, economic evaluation, plant location, and 
site preparation and structures. The volume is arranged 
so the more detailed material, such as plant layout and 
location, may be omitted without hindering the presenta- 
tion of other subject matter. 

All cost data have been brought up-to-date, and quite 
a few references dealing with recent developments in cost 


estimation techniques and equipment selectiun have been 
added since the last edition. All references for collateral 
reading have been placed together at the end of the text 
material, rather than at the end of each chapter. 

The rearrangement of the order of presentation of 
material undoubtedly has proxided a substantial improve- 
ment in the book ; however, there are a few specific 
criticisms, Of the topics presented, the section on economic 
evaluation is least satisfying, since the general subject is 
dealt with at considerable length but only the most simple 
concepts are presented quantitatively. A description of the 
economics of selecting alternates would be of particular 
value, as would be a more detailed consideration of the 
time value of money in economic calculations. The one 
entirely new section included in the fourth edition deals 
with nuclear chemical plant design, covering primarily 
shielding and criticality calculations. Both subjects are 
treated more quantitatively and with more detail elsewhere. 

The reviewer feels that the present edition is much 
clearer in presentation of material than previous editions. 
the thoroughness in coverage of cost estimation methods 
and data, together with the logical and clear presentation 
of the general problem of chemical plant design, make the 
book one to be recommended. J. B. Burr 
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